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Abstract: We study small hairy black holes in a consistent truncation of A/" = 8 gauged 

supergravity that consists of a single charged scalar field interacting with the metric and a 

C^h' U{1) gauge field. Small very near extremal RNAdS black holes in this system are unstable 

f— I ■ to decay by superradiant emission. The end point of this instability is a small hairy 

black hole that we construct analytically in a perturbative expansion in its charge. Unlike 

their RNAdS counterparts, hairy black hole solutions exist all the way down to the BPS 

^^ , bound, demonstrating that AA = 4 Yang Mills theory has an 0{N'^) entropy at all energies 

00 '. above supersymmetry. At the BPS bound these black holes reduce to previously discussed 

regular, supersymmetric horizon free solitons. We use numerical methods to continue the 

construction of these solitons to large charges and find that the line of soliton solutions 



^^ ' terminates at a singular solution 5 at a finite charge. We conjecture that a one parameter 

^-p ■ family of singular supersymmetric solutions, which emerges out of S, constitutes the BPS 

limit of hairy black holes at larger values of the charge. We analytically determine the near 
singularity behaviour of S, demonstrate that both the regular and singular solutions exhibit 
an infinite set of damped 'self similar' oscillations around S and analytically compute the 
C^ I frequency of these oscillations. At leading order in their charge, the thermodynamics of 

the small hairy black holes constructed in this paper turns out to be correctly reproduced 
by modeling these objects as a non interacting mix of an RNAdS black hole and the 
supersymmetric soliton in thermal equilibrium. Assuming that a similar non interacting 
model continues to apply upon turning on angular momentum, we also predict a rich 
family of rotating hairy black holes, including new hairy supersymmetric black holes. This 
analysis suggests interesting structure for the space of (yet to be constructed) hairy charged 
rotating black holes in AdS^ x S^, particularly in the near BPS limit. 
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1. Introduction 

Black hole solutions of JIB theory on AdS^ x S^ constitute the thermodynamic saddle 
points of AA = 4 Yang Mills theory on S^ via the AdS/CFT correspondence. A complete 
understanding of the space of stationary black hole solutions in AdS^ x 5^ is consequently 
essential for a satisfactory understanding of the state space of AA = 4 Yang Mills theory 
at energies of order A^^. While the Kerr RNAdS black hole solutions are well known 
in, ^, ^, ^, ^, ^, 0, it seems likely that several additional yet to be determined families of 
black hole solutions will play an important role in the dynamics and thermodynamics of 
TV" = 4 Yang Mills theory. 

In this paper we will construct a new class of asymptotically AdS^ x S^ black hole 
solutions. The black holes we construct are small, charged, and have parametrically low 
temperatures; our construction is perturbative in the black hole charge. Our solutions 
are hairy, in the sense that they include condensates of charged scalar fields^. In the 
BPS limit these hairy black holes reduce to regular horizon free solitons. We also use 
numerical techniques to continue our perturbative construction of these solitons to charges 
of order unity, and uncover an intricate self similar behaviour in the space of solitons in 
the neighborhood of a finite critical value of the charge. The solutions presented in this 
paper suggest a qualitatively new picture of the near BPS spectrum of A/" = 4 Yang Mills. 

Our perturbative construction of hairy black holes is close in spirit and technique to the 
constructions presented in the recent paper |l^], which may be regarded as an immediate 
precursor to the current work. For this reason we first present a brief review of |14] before 




, ]lu, hi, |l2| and references therein for reviews of recent work - sparked by an observation by 
on hairy black branes in AdS spaces. 



turning to a description and discussions of the new black hole solutions constructed in this 
paper. 



It was demonstrated in |14| that small charged black holes in global AdS spaces are 
sometimes unstable to the condensation of charged matter fields. More precisely any system 
governed by the Lagrangian 
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possesses small RNAdS black holes that are unstable to decay by superradiant discharge of 
the scalar field (/> whenever e > A. The end point of this superradiant tachyon condensation 
process is a hairy black hole. The authors of [14| constructed these hairy black hole solutions 
(working with a particular toy model Lagrangian of the form ( |1.1[ )) in a perturbation 
expansion in their mass and charge. At leading order in this expansion, the hairy black 



holes of 1 14 1 are well approximated by a non interacting mix of a small RNAdS black hole 



and a weak static solitonic scalar condensate. In particular, it was shown in |14| that the 
leading order thermodynamics of small hairy black holes could be reproduced simply by 
modeling them as a non interacting mix of an RNAdS black hole and a regular charged 
scalar soliton. 



The results of [14| suggest that the density of states of certain field theories with a 
gravity dual description might be dominated in certain regimes by previously unexplored 
phases consisting of an approximately non interacting mix of a normal charged phase and 
a Bose condensate. In order to make definitive statements about the actual behaviour of 
AA = 4 Yang Mills theory, however, it is necessary to perform the relevant calculations 
in IIB supergravity on AdS^ x S^ rather than a simple toy model Lagrangian; this is the 
subject of the current paper. As the IIB theory on AdS^ x 5"^ is a very special system, the 
reader might anticipate that hairy black holes in this theory have some distinctive special 
properties not shared by equivalent objects in the toy model studied in [U] at least at 
generic values of parameters. As we will see below, this indeed turns out to be the case. 

In order to avoid having to deal with the full complexity of IIB SUGRA, in this paper 
we identify^ work with a consistent truncation of gauged J\f = 8 supergravity (itself a 
consistent truncation of IIB SUGRA on AdS^ x S^). Of the complicated spectrum of 
M = 8 supergravity, our truncation retains only a single charged scalar field </>, a gauge 



^In unpublished work, S. Gubser, C. Herzog and S. Pufu have independently identified this consistent 
truncation, and have numerically investigated hairy black branes in this set up. We thank C. Herzog for 
informing us of this. 



field A^ and the metric.^ The Lagrangian for our system is given by 

' 8 '^ 16 V 4 (4 + (/.</.*) 
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and is of the general structure ( |1 . 1| ) with e = 2 and A = 2. As e = A, small very near 
extremal black holes in this system lie at the precipice of the super radiant instability 



discussed in |14|. The question of whether they fall over this precipice requires a detailed 
calculation. We perform the necessary computation in this paper. Our results demonstrate 
that very near extremal small black holes do suffer from the super radiant instability; we 
proceed to construct the hairy black hole that constitutes the end point of this instability. 

The calculations presented in this paper employ techniques that are similar to those 
used in [^. In particular we work in a perturbation expansion in the scalar amplitude 
using very near extremal vacuum RNAdS black holes as the starting point of our expansion. 
This expansion is justified by the smallness of the charge of our solutions. We implement 
our perturbative procedure by matching solutions in a near (horizon) range, intermediate 
range and far field range. This matching procedure is justified by the parametrically 
large separation of scales between the horizon radius of the black holes we construct and 
AdS curvature radius. We refer the reader to the introduction of [0] for a more detailed 
explanation of physical motivation for this perturbative expansion and its formal structure. 
Actual details of our calculations may be found in sections |3| and |^ below. In the rest of 
this introduction we simply present our results and comment on their significance. 

We study black holes of charge Q ^ (normalized so that each of the three complex Yang 

Mills scalars has unit charge) and mass M (normalized to match the scaling dimensions of 

dual operators). We find it convenient to deal with the 'intensive' mass and charge, m and 

q, given by 

M Q 



As in 1 14 1, in this paper we are primarily interested in small black holes for which q <^1 
and m <^1. We now recall some facts about RNAdS black holes in this system. First, the 
masses of such black holes obey the inequality 

Black holes that saturate this inequality are extremal, regular and have finite entropy (see 



subsection 6J for more details). The chemical potential HextiQ) of extremal black holes is 

given by 

fiextiq) = l + 2q-6q^ + 0{q^) 

^Under the AdS/CFT correspondence, (j)_ is dual to the operator TrX"^ +TrY'^ +TrZ'^ while A^ is dual 
to the conserved current J^^ + J^^ + J^^ ■ Here X, Y and Z denote the three complex chiral scalars in 
the N = A Lagrangian. 

"^Q is the charge of the black hole solutions under each of the three diagonal U(l) Cartan's of SO(6). 
The consistent truncation of this paper forces these three U(l) charges to be equal. 



and approaches unity in the hmit of smah charge. Note also that extremal black holes lie 
above the BPS bound 

rriBPsiq) = 3g 

and in particular that mextio.) — "mBPsil) = 3g^ + 0{q^) 

We will now investigate potential superradiant instabilities (see the introduction of 



1 14] for an explanation of this term) of these black holes. Recall that a mode of charge e 



and energy a; scatters off a black hole of chemical potential /x in a superradiant manner 



whenever fxe > co. The various modes of the scalar field (f) in (1^) have energies w = 2, 3, . . . 
and all carry charge e = 2. As the chemical potential of a small near extremal black hole 
is approximately unity, it follows that only the ground state of 4> (with uj = A = 2) could 
possibly scatter of a near extremal RNAdS black hole in a superradiant manner. This mode 
barely satisfies the condition for superradiant scattering; as a consequence we will show 
in this paper that small RNAdS black holes in (|l]^) suffer from a superradiant instability 
into this ground state mode only very near to extremality, i.e. when 

m-mext{q)<Qq^ + 0{q^). 

Unstable black holes eventually settle down into a new branch of stable hairy black hole 
solutions. We have constructed these hairy black holes in a perturbative expansion in their 
charge in Section |^; we now proceed to present a qualitative description of these solutions 
and their thermodynamics. 

Recall that the zero mode of the scalar field (/> obeys the BPS bound (and so is super- 
symmetric) at linear order in an expansion about global AdS^. It has been demonstrated 
in jl^, |l^, |l^ (and we reconfirm in section ^ below) that this linearized BPS solution 
continues into a nonlinear BPS solution upon increasing its amplitude. In this paper we 
will refer to this regular solution as the supersymmetric soliton. The hairy black holes of 
this paper may approximately be thought of as a small, very near extremal RNAdS black 
hole located in the center of one of these solitons. Although the soliton is supersymmetric, 
the black hole at its center is not, and so hairy black holes are not BPS in general. These 
solutions exist in the mass range 

3q<m<3q + 3q'^ + 0{q'^) (1.3) 

At the lower bound of this range ( |7.14| ) hairy black holes reduce to the supersymmetric 
soliton. At the upper bound (which is also the instability curve for RNAdS black holes) 
they reduce to RNAdS black holes. 

In Fig. |l] below we have plotted the near extremal micro canonical 'phase diagram' for 
our system. As is apparent from Fig. Q our system undergoes a phase transition from an 
RNAdS phase to a hairy black hole phase upon lowering the energy at fixed charge. This 
phase transition occurs at the upper end of the range (|l.3| ). Note that the phase diagram 
of Fig. H has several similarities with the phase diagram depicted in Fig. 1 in |1J] ; however 
there is also one important difference. The temperature of the hairy black holes of this 
paper decreases with decreasing mass at fixed charge, and reaches the value zero at the 




Figure 1: Phase diagram as a function of charge q (x axis) and mass m (y axis) at small q. The 
solid blue line at the bottom is the BPS bound along which the soliton lives. Hairy black holes 
exist - and are the dominant phase - in the shaded region. RNAdS black holes are the only known 
solutions (so in particular the dominant phase) in the unshaded region above the solid red curve at 
the top. RNAdS black holes also exist (but are dynamically unstable and thermodynamically sub 
dominant) between the solid red curve and the dashed curve. The solid red curve is described by 
ra — ?>q-\- 3g^ + C'(g'*), while the blue curve hy m — 3g. The dashed curve corresponds to extremal 
RNAdS black holes and is given by m = 3g + iq^ — Qq^ + 0{q'^). The curves have not been drawn 
to scale to make the diagram more readable. 



BPS bound. In contrast the temperature of the hairy black holes of [14| increases with 
decreasing mass (at fixed charge), approaching infinity in the vicinity of the lower bound. 
As we have emphasized, the phase diagram depicted in Fig. |l| apphes only in the limit 
of small charges and masses. We would now like to inquire as to how this phase diagram 
continues to large charges and masses. In order to address this question we first focus 
on solitonic solutions. These solutions may be determined much more simply than the 
generic hairy solution, as they obey the constraints of supersymmetry rather than simply 
the equations of motion. It turns out that spherically symmetric super symmetric solutions 
are given as solutions to a single nonlinear, second order ordinary differential equation 



1 15, 1£, 17]. The solitons constitute the unique one parameter set of regular solutions 
to this equation. It is easy to continue our perturbative construction of the solitonic 
solutions to large charges by solving this equation numerically: in fact this exercise was 



already carried out in |16|. This numerical solution reveals that the solitonic branch of 
solutions terminates at a finite charge qc = 0.2613. For q > Qc there are no supersymmetric 
spherically symmetric solutions to the equations of motion of ( |1.2| )^. 

Recall that solitons constitute the lower edge of the space of hairy black hole solutions 
of Fig. ||. The non existence of regular supersymmetric solutions for q > Qc might, at first, 
suggest that at these charges the space of hairy black hole solutions terminates at a mass 
greater than 3q (i.e. does not extend all the way down to supersymmetry). While this is 
a logical possibility, we think it is likely that the truth lies elsewhere. As we will explain 



^To be more precise, there are smooth solitonic solutions up to a slightly higher value qm = 0.2643, but 
in a sense that will be explained in section H, the point qc marks the boundary between regular solitonic 
solutions and singular ones. 



in section ^ the solitonic branch of supersymmetric solutions terminates in a distinguished 
singular solution S. It turns out that S is also the end point (or origin) of a one parameter 
set of supersymmetric solutions that are all singular at the origin. The charges of these 
solutions increase without bound (indeed we have found an explicit analytic solution for 
the singular supersymmetric solution in the limit of arbitrarily large charge) . The two one 
parameter families of solutions, regular and singular ones, are joined at the special solution 
S. We conjecture that smooth hairy black hole solutions exist in our system at every q 
and for m > 3q. Upon taking the limit ttt, — )■ 3g, these smooth solutions reduce to the 
smooth soliton for q < Qc but reduce to the singular supersymmetric solutions described 
above when q> qc. In summary, we conjecture that the phase diagram of our system takes 
the form displayed in Fig. below. 




Figure 2: Conjectured phase diagram as a function of charge q (x axis) and mass m (y axis) for 
all values of q. The blue line at the bottom is the BPS bound along which the regular soliton lives 
(straight part) and the singular supersymmetric solutions (wiggly part). The solid red curve at the 
top marks the phase transition between the regime of RNAdS black holes (above the line) and that 
of smooth hairy black holes (below). The black curve indicates a phase transition between different 
types of hairy black holes. This curve originates from the BPS line at the black dot which is close 
to the point q = Qc and could end either in the bulk of the hairy black hole region or could extend 
all the way up to the red line. 

The distinguished solution S clearly plays a special role in the space of spherically 
symmetric supersymmetric solutions. In subsection |5.2| we analytically determine the near 
singularity behaviour of this solution. Viewing the 2nd order differential equation that 
determines supersymmetric solutions as a dynamical system in the 'time' variable Inr, we 
demonstrate that the solution 5 is a stable fixed point of this system, and analytically 
compute the eigenvalues that characterize the approach to this fixed point. This eigen- 
value has an imaginary part (which damps fluctuations) and a real part (that results in 
oscillations). Solitonic - and singular - solutions in the neighborhood of S may be thought 
of as configurations that that flow to S at large Inr. The oscillations^ referred to above 
result in the following phenomenon: the system develops a multiplicity of supersymmetric 



®We are extremely grateful to M. Rangamani for suggesting that we look for this 'self similar' structure 
in the space of solitons in the neighborhood oi q = qc- The results reported in this paragraph are the 
outcome of investigations that were spurred directly by this suggestion. 



solitonic (or singular) at charges q when q comes near enough to qc- The number of solu- 
tions diverges as q ^ qc- The space of solitonic and singular supersymmetric solutions are 
usefully plotted as a curve on a plane parametrized by the charge q and the expectation 
value of the operator dual to the scalar (p. On this plane supersymmetric solitons spiral 
into the point S, while the singular solutions spiral out of the same point (see Fig. 15); 
the two spirals are non self intersecting ^. We find this extremely intricate structure of 
supersymmetric solutions quite fascinating, and feel that its implications for AA = 4 Yang 
Mills physics certainly merits further investigation. 

In this paper we have so far only considered charged black holes with vanishing angular 
momentum. Such solutions are spherically symmetric; i.e. they preserve the 50(4) = 
SU{2) X SU{2) rotational isometry group. As all known supersymmetric black holes in 
AdS^ X S^ possess angular momentum ||l^, |l^, |2^], it is of interest to generalize the study 
of this paper to black holes with angular momentum. Let us first consider a spinning 
Kerr RNAdS black hole that preserves only [/(I) x C/(l) e SU{2) x SU{2). Perturbations 
about such a solution are functions of an angle and a radius and are given by solutions 
to partial rather than ordinary differential equations. However there exist RNAdS black 
holes with self dual angular momentum. The angular momentum of such a black hole lie 
entirely within one of the two SU{2) above, and so preserve a U{1) x SU{2) subgroup of the 
rotation group. Perturbations around these solutions may be organized in representations 
of SU{2) and so obey ordinary rather than partial differential equations.^ Consequently, 
the generalization of the hairy black hole solutions determined in this paper to solutions 
with self dual spin appears to be a plausibly tractable project; which we however leave to 
future work. 

Even though we do not embark on a serious analysis of charged rotating black holes 
in this paper, we do present a speculative appetizer for this problem. In more detail we 
present a guess (or a prediction) for the leading order thermodynamics of these spinning 
hairy solutions. Our guess is based on the observation that the thermodynamics of the hairy 
solutions constructed in this paper can be reproduced, at leading order, by modeling the 
hairy black hole as a non interacting mix of a RNAdS black hole and the soliton. In section ^ 
we simply assume that a similar model works for charged spinning black holes, and use this 
model to compute the thermodynamics of a certain class of spinning hairy black holes. The 
most interesting aspect of our results concern the BPS limit. Our non interacting model 
predicts that extremal hairy black holes are BPS at every value of the angular momentum 
and charge. This is in stark contrast with Kerr RNAdS black holes that are BPS only 
on a CO dimension one surface of the space of extremal black holes. According to our non 
interacting model, BPS hairy black holes are a non interacting mix of Gutowski Reall black 
holes [18, ^, 2C] and supersymmetric solitons. Such a mix is thermally equlibrated at all 
values of charge and angular momentum because of an important property of Gutowski 
Reall black holes; their chemical potential is exactly unity^. We find this result both 
intriguing and puzzling (see e.g. |21]). We emphasize that our prediction is based on 



^We are very grateful to V. Hubeny for suggesting the possibility of a spiral structure for these curves. 
*We thank A. Strominger for pointing this out to us. 
^We are very grateful to S. Kim for explaining this to us. 



the non interacting superposition model, which may or may apply to actual black hole 
solutions. We leave further investigation of this extremely interesting issue to future work. 

This paper has been devoted largely to the study of small very near extremal charged 
black holes in AdS^ that are smeared over S^. As uncharged small smeared black holes are 



well known to suffer from Gregory-Laflamme type instabilities |22|, the reader may wonder 
whether the black holes studied in this paper might suffer from similar instabilities. We 
believe that this is not the case. Recall that the likely end point of a Gregory Laflamme 
type instability is a small black hole of proper horizon radius r^ localized on the S^. In 
order that this black hole be near extremal, it has to zip around the S^ at near the speed of 
light, i.e. at V = 1 — 5 with (5 <C 1. The AdS^ charge of such a black hole is given by g oc -ft 
while its energy above extremality of such a black hole is given by m — 3q ex rJj\/6.^^ We 
may now solve for rn and 6 as a function of q and m — 3q. In the near BPS limit of interest 
to this paper m — ^q ^ q^ and we find rj^ oc (72 and \f5 oc ^Jq. It follows that the entropy 
of such a localized black hole oc r|^ oc (72^ 7 ^ and so is smaller than the entropy, (oc q^), 
of the black holes studied in this paper. As the the black holes studied in this paper have 
higher entropy than S^ localized black holes with the same charge, there seems no reason 
to expect them to suffer from Gregory- Laflamme type instabilities^^. Another pointer to 



the same conclusion is the fact it was very important for the analysis of |22] that the black 
holes they studied had negative speciflc heat. The charged black holes at the center of the 
hairy solutions here all have positive specific heat ^^. 

We also note that the Gubser Mitra instability [^, ^] afflicts three equal charge black 
holes only when the black holes in question have large enough charge. It follows that the 
small black holes primarily studied in this paper do not suffer from Gubser Mitra type 
instabilities^^. 

It is conceivable that the solutions presented in this paper might suffer from further 
superradiant instabilities, once embedded in IIB theory on AdS^ x S^ . In order to see 
why this might be the case, let us recall once again why the field (p - dual to the chiral 
Yang Mills operator TrX^ + TrY^ + TrZ'^ condensed in the presence of very near extremal 
charged RNAdS black hole. The reason is simply that the energy A = 2 of this field is 
equal to its charge e = 2. As a consequence the Boltzmann suppression factor, e~^^ ' 
of this mode exceeds unity when ^ > 1 causing this mode to Bose condense. However 
exactly the same reasoning applies to, for instance, the field (pn dual to the chiral operator 



""^"To see this let the sphere be given by equations |zip + |2;2p + jzap = 1 where Zi are the three complex 
embedding coordinates. A black hole we study is located at jzij = \z2\ = \z3\ = -js, and moves with speed 
^ ~^' in each of the three orthogonal planes. Let the proper mass of the black hole be m,p oc rjj. Its angular 

momentum in each plane, q — r x p, is given by ^ x '"''^V- — - = IHeli 1 where 7 = "*" The 

energy of the black hole is 771^7. 

^Were we interested in black holes with small q and m — 3g ~ 0(9) then we would have found 5 ~ 1 and 
rjf ~ q. The entropy of the localized black hole would then have been oc g 7 > g 2 . As such black holes can 
increase their entropy by condensing, they presumably do suffer from Gregory Laflamme type instabilities. 

^^We thank V. Hubeny and M. Rangamani for discussions on this point. 

'^^We thank M. Rangamani for a discussion on this point. 



TrX^ + TrY^ + TrZ^ all of which have A = e ^^. It seems likely that there exist other 
hairy solutions in which some linear combination of (pn (rather than simply (P2) condense 
^^ . It is important to know whether any solution of this form has higher entropy than the 
black holes with pure <j)2 condensate presented in this paper. If this is the case then the 
hairy black holes of our paper would likely suffer from superradiant instabilities towards the 
condensation to the entropically dominant black hole. On the other hand the black holes of 
this paper, with 4>2, the lightest chiral scalar operator that preserves all discrete symmetries 
of the problem, as the only condensate, are quite special. It seems quite plausible to us that 
the solution presented in this paper has the largest entropy of all the hairy solutions with 
Pn condensates. If this is indeed the case then the hairy black hole solutions presented in 
this paper constitute the thermodynamically dominant saddle point of AA = 4 Yang Mills 
very near to super symmetry; and the entropy of AA = 4 Yang Mills very near to the BPS 



bound is given the formula (6.17) below. 

To end this introduction we would like to emphasize that the black hole solutions of 
this paper give a qualitatively different picture of the density of states of AA = 4 Yang Mills 
theory at finite charge compared to a picture suggested by RNAdS black holes. As we have 
seen above, there exist no RNAdS black holes with masses between mBPsil) and mext{Q), 
a fact had previously been taken to suggest that, for some mysterious reason, there are less 
than 0{N'^) states in Yang Mills theory between rriBPsio) and m^xtiQ)- The new black 
hole solutions of this paper establish, on the other hand, that AA = 4 Yang Mills theory has 
0{N'^) states all the way down to the BPS bound at least at small charge, and plausibly at 
all values of the charge (see Fig. |2|).^^ The saddle point that governs near BPS behaviour is 
a mix of a charged Bose condensate and a normal charged fluid. It would be fascinating to 
find some (even qualitative) confirmation of this picture from a direct field theory analysis. 

2. A Consistent Truncation and its Equations of Motion 

2.1 A Consistent Truncation of Gauged Supergravity 

A/" = 8 gauged supergravity constitutes a consistent truncation of IIB theory on AdS^ x S^ . 
In addition to the metric, the bosonic spectrum of this theory consists of 42 scalar fields, 
15 gauge fields and 12 two form fields. The scalars transform in the 20 +10c +1 +1 of 
S0(6), the gauge fields transform in the 15 dimensional adjoint representation, while the 
two form fields transform in the 6c representation of 50(6). 

It has been shown |^] that AA = 8 gauged supergravity admits a further consistent 
truncation that retains only the scalars in the 20 and the vector fields in the 15 together 



^^This statement is true more generally of every operator in the A/" = 1 chiral ring of the theory. 

^^In the BPS limit any linear combination of pnS can condense and we have an infinite dimensional moduli 
space of solutions (see pa, [L7|). We expect the introduction of a black hole to lift this moduli space, to a 
discrete set of solutions. 

^®This difference is starkest in the limit of large charge, i.e. in the Poincare Patch limit. The energy 
density, pE, of RNAdS black branes is bounded from below by cpl where pQ is the charge density. Fig. H 
, on the other hand predicts that the energy density of a charged black brane can be arbitrarily small at 
any given value of the charge density. 
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with the metric, setting all other fields to zero. The action for this consistent truncation 
is given by 



1 / 2 ^^ ^'-'^ 

where 

(2-2) 

Here (i, j, • • • ) denote the 5*0(6) vector indices and {fi,i^,- • •) are the space time indices. 
Tij are symmetric unimodular (i.e. Tij is a matrix of unit determinant) SO{6) tensors. 
Further A^ is the rank of the gauge group of the dual AA = 4 Yang Mills theory, and we 
work in units in which the AdS^ with unit radius solves ( p.l[ ). 

We will now describe a further consistent truncation of (p^) • For this purpose we find 
it useful to move to a complex basis for the 50(6) vector indices that appear summed in 



(2.1). Let {xj j = !,••• ,6) denote 50(6) Cartesian directions. We define the complex 



coordinates 

X2J-1 + iX2j = Zj, X2J-1 - iX2j = Zj j = 1, 

We will now argue that the restriction 



-Lzizi — J^Z2Z2 — -'2323 

V4 + <j)(t)* (2-3) 



T- - — T- - — T- - — —- 

-l-ZiZl — ^Z2Z2 — -^2323 ~ A 



-t^izi — ''-Z2Z2 — -'-z^z-i 



AZlZl _ \Z2Z2 _ AZ3Z3 _ r)- A 



All Others =0 

constitutes a consistent truncation of (p.l|). To see this is the case note that the permu- 
tations of labels (1,2,3), as also separate rotations by vr in the zi, Z2 and 23 planes, can 
each be generated by separate 5*0 (6) gauge transformations. It follows that these discrete 



transformations are symmetries of (p.l|). Now it is easy to convince oneself that (2^) is 
the most general field configuration of (|2.1| ) that is invariant separately under each of these 
four discrete symmetries. It follows that ( |2.3| ) is a consistent truncation of the system (p.l|). 
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The consistent truncation (|2.3| ) is governed by the Lagrangian 
1 



iV2 



V^ 



A^ 1^ 



2 ^ ^ 8 ^ 16 V 4 (4 + (/.(/.*; 



-{R+12)--F,^F^'^--i\D 



16 



4(4 + 00*) 



■400' 



-0^0 = 9^0 - 2M^0 



(2.4) 



Note that has charge 2 and m^ = —4. Under the AdS/CFT dictionary this field maps 
to an operator of dimension A = 2. Note also that the kinetic term of the gauge field has 
the factor prefactor | rather than the (more usual) 4 as employed, for instance, in |14]. ^^ 

2.2 Equations of Motion 

We now list the equations of motion that follow from varying (2.4). We find the Einstein 
equation 



Rfiy - i^a^iuR - Qgiiu - -i;T^u + o^/ 



nmat 



(2.5) 



where 



rpEM _p ap _± Tpaap 



TP"' =- [D^cp {D,cl>r + Z?.0 (D^cPT] - -g^,AD^ct>\^ + 200' 



a^lv 



(2.6) 



4(4 



9^(00*)9,(00*) - -<?;..[a,(00*)]^ 



the Maxwell equation 



V^F^- = ^ [<t>{D,^ 



*D, 



and the scalar equation 



*M2 



4(4 + 00*)2 



V= 



2(4 



+ 4 



= 0. 



(2.7) 



(2.8) 



'^Consequently gauge fields and chemical potentials in this paper and ||l4[| are related by 



^/lere \/ ^^there: l-^'here 



ij f^there 



Note also that Gs was set to to unity in |14|, while Gs ~ ySr in this paper. It follows that 



Mhe 



iV2 



_^ i^r Ohere _£ „ 

TV N-' TV 



^ here ^ / /I 

Try ^ \ Tl^there 



The factor of J^ above ensures that TdSthere = dMthere — jJ-theredQthere implies TdShere = dMhere 

Sfj,heredQhere as is required on on physical grounds. 
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In this paper we study static spherically symmetric configurations of the system (| 
We adopt a Schwarzschild like gauge and set 

ds'^ = -f{r)dt^ + g{r)dr'^ + r'^d^l 
At = A{r) 

Ar = ^,- = 



(2.9) 



The four unknown functions /(r), g{r\ A{r^ and (/)(r) are constrained by Einstein's 
equations, the Maxwell equations and the scalar equations. It is possible to demonstrate 
that /, 5, A^ (p are solutions to the equations of motion if and only if 



E^ 



Eo 




(2.10) 



^3 =2A{r)g{r)(l){rf + A'{r) 



f'jr) ^ g'{r) 6 
f{r) g{r) r 



Ea 



4 + . 



where 



V^<h(r) 



V^(h(r) 



9[r) 



1 + 



A{7 



- 2A"{r) 
cP'irf 



m 



g[r) 



{rf+4Y 



f'jr) 
fir) 



+ 



(t)'{r) + 2(t)"{r) -g'{r)(l)'{r 



2g{ry 

The equations Ei and E2 are derived from the rr and tt components of the Einstein 
equations, E^, is the t component of the Maxwell equation and E/^ is the equation of the 
scalar field. 



As in [14| the equations (2.10) contain only first derivatives of / and g, but depend on 
derivatives up to the second order for and A. It follows that ( |2.10| ) admit a 6 parameter 
set of solutions. One of these solutions is empty AdSs space, given by /(r) = r^ + 1, 
S'l'^) = T+r^' ^('") ~ *^('") ~ ^- ^^ ^"^^ interested in those solutions to ( 2.10D that asymptote 
to AdS space time, i.e. solutions whose large r behaviour is given by 

/(r) = r2 + l + 0(l/r2) 
1 



9{r} 



+ 0(l/r^) 



l + r2 

A{r) = 0{l) + 0{l/r^) 
0(r) = 0(l/r2) 



(2.11) 
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As in ||T^ it turns out that these conditions effectively impose two conditions on the solu- 
tions of ( |2.1C1| ) , so that the system of equations admits a four parameter set of asymptot- 
ically AdS solutions. We usually also be interested only in solutions that are regular (in 
a suitable sense) in the interior. This requirement will usually cut down solution space to 
distinct classes of two parameter space of solutions; the parameters may be thought of as 
the mass and charge of the solutions. 

2.3 RNAdS Black Holes 

The AdS-Reissner-Nordstrom black holes constitute a very well known two parameter set 
of solutions to the equations ( ^.10 ). These solutions are given by 



fir) =^ - ^ ^2 + ^ + 1 

=^Ar'-R'][r' + r\R' + i)-f^'R'] 

q(r) = J_ (2.12) 

9[r) ^(^^ 

A{r) =/x (l - ^ 
6(r) = 



where fi is the chemical potential of the RNAdS black hole. The function V{r) in ( 2.12| ) 



vanishes at r = i? and consequently this solution has a horizon at r = R. In fact, it can 
be shown that R is the outer event horizon provided 



fi'' <{1 + 2R^). (2.13) 



As explained in [p!^ ] and in the introduction, ( |2.12 ) is unstable to superradiant decay 



provided in the presence of field of charge e and minimum energy A provided e// > A. Now 
our field (p has A = 2 and e = 2. Moreover, in the limit i? — )• 0, RNAdS black holes have 
fj, < 1 (this inequality is saturated at extremality) . It follows that small extremal black 
holes lie at the edge of instability, as mentioned in the introduction. We show below that 
very near extremal RNAdS black holes do in fact suffer from super radiant instabilities. 

3. The Supersymmetric Soliton in Perturbation Theory 

In this section we will construct the analogue of the ground state soliton in |14]. The 



new feature in here is that the soliton turns out to be supersymmetric (this is obvious at 
linearized order). 

In this section we generate the solitonic solution in perturbation theory. We use only 
the equations of motion without imposing the constraints of supersymmetry, but check 
that our final solution is supersymmetric (by verifying the BPS bound order by order in 
perturbation theory). This method has the advantage that it generalizes in a straightfor- 
ward manner to the construction of non supersymmetric hairy black holes in subsequent 
sections. 
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In Section ^ we will revisit this solitonic solution; we will rederive it by imposing the 
constraints of supersymmetry from the start. That method has the advantage that it 
permits a relatively simple extrapolation of super symmetric solutions to large charge. 

3.1 Setting up the perturbative expansion 

We now turn to the description of our perturbative construction. To initiate the perturba- 
tive construction of the supersymmetric soliton we set 

n 
^ oo 

9{r) = -—-,+Y.''''92n{r) 
1 + r^ ^-^ 

ra=l 

oo 

yl(r) = l + ^62M2„( 

n=l 

oo 



ra=l 

(3.1) 

n=l 



n=l 



and plug these expansions into ( 2.10| ). We then expand out and solve these equations order 
by order in e. All equations are automatically solved up to 0{e). At order e^" the last 
equation in ( |2.1C1| ) is trivial while the first three take the form 

^(r=(l+r=)V2„w)=P£'(r) 

^^^■^-^^i±^...« + P.T« (3.2) 



dr \1 + r^ J r 

d f , dA2nir) \ _ (A) 



On the other hand, at order 6^"+^ the first three equations in ( 2.10| ) is trivial while the last 
equation reduces to 

iilT^Tr [d + -')^-+i(-)] ) = P^tUr) (3.3) 

Here the source terms P2n('^)' ^2ni''")^ ^2n i''^) ^^^'^ P2„+i(^) ^^^ ^o^P^^t^^y determined 
by the solution to lower orders in perturbation theory, and so should be thought of as 
known functions, in terms of which we wish to determine the unknowns /2n, g2n, ■^2n and 

4>2n+l- 

The equations ( |3.2| ) are all easily integrated. It also turns out that all the integration 
constants in these equations are uniquely determined by the requirements of regularity. 
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normalisability and our definition of e, exactly as in |14]. The solution is given by 

1 



92n{r) 



Pii{s)ds 



2(1 + 2s^ 



/2n(r)=-(l+r2) / 

J r 

92n^\\r) = -- I as 



92n{s) + P^iHs) 



M 



ds 



(3.4) 



Prnilis')ds' 



3.2 The Soliton up to 0{e^) 

The perturbative procedure outlined in this subsection is very easily implemented to arbi- 
trary order in perturbation theory. In fact, by automating the procedure described above, 
we have implemented this perturbative series to 17th order in a Mathematica programme. 
In the rest of this subsection we content ourselves with a presentation of our results to 



92[r) 
f2{r) 

A2{r) 

54 (r) = 

h{r) = 

A,{r) = 

Ur) = 







4(1 + r2) 


1 


8(1 + r2) 


1 


8(l+r2)3 


^4 


192(1 +r2)5 


^4 


192(1 + r2)5 


r^ + 3r2 + 3 


384(r2 + l)3 


g^4 _^ 4^2 _^ 55 



(3.5) 



(3.6) 



2304 (r2 + l) 



Aeir) 

(j)j{r) 



r^ {Ar^ + 15r2 + 20) 
7680 (r2 + l)'^ 
12r8 + 45r6 + 60r'' + 20r2 + 5 

23040 (r2 + 1)5 
6r^ + 30r^ + 60r^ + 55r^ + 25 
23040 (r2 + if 
120r^ + 460r6 + 1095r'' + 558r2 + 2368 
460800 (r2 + l)^ 



(3.7) 
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r4 



fsir) 



Asir) 



r- (I69r^ + 1024r^ + 2640r'^ + 3320r2 + 2180) 
2211840 (r2 + l)^ 
(5 (l69rS + 1024r*5 + 2580r^ + 3344r2 + 2288) r^ + 3096) r^ + 516 

11059200 (r2 + l)^ 
(5 ((169 (r^ + 7r2 + 2l) r^ + 5819) r^ + 5543) r^ + 14721) r^ + 4191 

22118400 (r2 + l)'^ 



(5 (I014r8 + 6124r6 + 18257r^ + 30484r2 + 36676) r^ + 75784) r^ + 155759 
(r) = -^^ — ^^ 



132710400 (r2 + 1)^ 



(3.8) 



The soliton obeys the BPS relation m = 3^ to the order to which we have carried out 
our computation (we present more details of the thermodynamics in Section |6|). 



4. The Hairy Black Hole in Perturbation Theory 

4.1 Basic Perturbative strategy 

We will now present our perturbative construction of hairy black hole solutions. In order 
to set up the perturbative expansion we expand the metric gauge field and the scalar fields 
as 



/(0=E^'"/2n(r) 

n=0 

oo 

n=0 

oo 

^(0=E^'"^2n(r) 

n=0 

oo 

:5;62"+V2n+l(r) 



(4.1) 



n=0 



where the unperturbed solution is taken to be the RNAdS black hole 
Ur,R)=V{r), go{r,R) ^ 



V{r) 

Ao{r,R) = ,io{l-^) (4.2) 

V{r) = l + r^(l- ^Vo + f + R' ^ ^ 
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The chemical potential of our final solution will be given by an expression of the form 

(4.3) 







n=0 




oo 




f^2n{R) ■■ 


= ^ ^J'{2n 
fc=0 


n2k 
,2fc)-K 


^(0,0) '■ 


= 1 





Note that, at the leading order in the perturbative expansion, /i = 1. 

Our basic strategy is to plug the expansion ( [4.ip into the equations of motion and 
then to recursively solve the later in a power series in e. We expand our equations in a 
power series in e. At each order in e we have a set of linear differential equations (see below 
for the explicit form of the equations), which we solve subject to the requirements of the 
normalisability of (/>(r) and /(r) at infinity together with the regularity of (/>(r) and the 
metric at the horizon. These four physical requirements turn out to automatically imply 
that A{r = R) = i.e. the gauge field vanishes at the horizon, as we would expect of a 
stationary solution. These four physical requirements determine 4 of the six integration 
constants in the differential equation, yielding a two parameter set of solutions. We fix 
the remaining two integration constants by adopting the following conventions to label 
our solutions: we require that (p{r) fall off at infinity like -% (definition of e) and that 
the horizon area of our solution is 27r^i?^ (definition of R). This procedure completely 
determines our solution as a function of R and e. We can then read of the value of /i in 
(|4.3| ) on our solution from the value of the gauge field at infinity. 

As in |]lj], the linear differential equations that arise in perturbation theory are difficult 
to solve exactly, but are easily solved in a power series expansion in R, by matching near 
field, intermediate field and far field solutions. At every order in e we thus have a solution 
as an expansion in R. Our final solutions are, then presented in a double power series 
expansion in e and R. 

4.2 Perturbation Theory at 0{e) 

In this section we present a detailed description of the implementation of our perturba- 
tive expansion at 0{e). The procedure described in this subsection applies, with minor 
modifications, to the perturbative construction at ©(e^™"*"^) for all m. 

Of course all equations are automatically obeyed at 0{e^). The only nontrivial equation 
at 0{e) is D'^4i = where D^ = V^ — 2iA^ is the linearized gauge covariantised derivative 
about the background ( |4.2| ). We will now solve this equation subject to the constraints of 
normalisability at infinity, regularity at the horizon, and the requirement that 

at large r. 
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4.2.1 Far Field Region (r > R) 

Let us first focus on tlie region r ^ R. In tiiis region the black hole ([4 .21) 



2 



/oM=^:4 ^2 — ^ — + ^ +1 

ffo(r) 



(4.4) 



m 

Ao{r) =/io (^1 - ^ 

oo 

lJ-0 =2_^R /^(0,2fc) 
fc=0 

/^(0,0) =1 

is a small perturbation about global AdS space. For this reason we expand 

oo 

<Arw = E^''C2fe)W' (4-5) 

fc=0 

where the superscript out emphasises that this expansion is good at large r. In the limit 
i? — )• 0, ( ^ ) reduces to global AdS space time with At = I. A stationary linearised 
fluctuation about this background is gauge equivalent to a linearised fluctuation with time 
dependence e~** about global AdS space with ^4^ = {At is the temporal component of the 
gauge field). The required solution is simply the ground state excitation of an m^ = —4 
minimally coupled scalar field about global AdS and is given by 

Co)(-) = T^ (4-6) 

The overall normalisation of the mode is set by our definition of e which implies 



We now plug ( [4.5[ ) into the equations of motion D (p = and expand to 0{R ) to solve 

out 
(1,2)- 

p2)o«t ^ (l)2)««t + R^{Dl)'^^ + ... 



for <t)"(i\\- Here D^ is the gauge covariant Laplacian about the background (|4.4|). Now 



where (Z?q)°"* is the gauge covariant Laplacian about global AdS space time with back- 
ground gauge field At = 1. It follows that, at 0{R'^), 

/■ 1-^2 \ out I out / j-\2\ out tout / j-\2\out 

[Do) (A(i,2) = -(^2) (f>{i,o) = -{D2) 



l + r2 



This equation is easily integrated and we find 

C2)(-) = - (yT^) (^ - 2 [^(0,2) - 2] log(r) + [/.(o,2) - 2] log (r^ + l) + -A_) 

(4.7) 
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We could iterate this process to generate (p'(i^2k) ^pto any desired order k. The equa- 



tions we need to solve, at order 0{R ), takes the form 



d 



d 



dr\l-\-r^dr 



[l + r' 



lOUt 

(l,2fc) 



(r) 



PiiU') 



(4.8) 



where P(T\j.\ {r) is a source function, whose form is determined by the results of the expan- 
sion at lower orders in perturbation theory. 

As in (^]^), it turns out that the expressions <;/>?" 2^1 ^"^^ increasingly singular as r — )• 
0. In fact it may be shown that the most singular piece of <P1i\i^\ scales like ^^, upto 
logarithmic corrections. In other words the expansion of </>°"* in powers of B? is really an 
expansion in -^ (upto log corrections) and breaks down at r ~ i?. 

To end this subsection we summarize our results to 0{E?'). We have 
1 



lOUt 



l + r- 



: + 



1 



l + r" 



A''(0,2) 



2 [;U(o,2) - 2] log(r) + [^(0,2) - 2] log (r^ + l) + 



+ 



R' 



H + 1 



(4.9) 



Expanding (p" (r) in a Taylor series about r = we find 



6°"*(r) 



[l-r2 + 0(r^)] 



+ R^ 



.2 
^(0,2) 



+ (^(0,2)-2)(21og(r) + l) + O(r2) 



(4.10) 



,R 



+ Oi-^){l + Oir')) + ... 



Note that this result depends on the as yet unknown parameter A*(o,2)- This quantity 
will be determined below by matching with the intermediate field solution of the next 
subsection. 

4.2.2 Intermediate Field Region r <C 1 and (r — i?) » R^ 

Let us now turn to intermediate region R^ <^ r — R <^ 1. Over these length scales the small 
black hole is far from a small perturbation about AdSs space. Instead the simplification 
in this region stems first from the fact that we focus on radial distances of order R ( 
r ~ i? ^ 1). Over these small length scales the background gauge field, which is of order 
unity, is negligible compared to the mass scale set by the horizon radius 4. 

A second simplification results from the fact that we insist that {r — R) ^ i?^, i.e. 
we do not let our length scales become too small. At these distances the back hole that 
we perturb around are effectively extremal (rather than slightly non extremal) at leading 
order. Moreover the black hole may also be thought of (at leading order) as a small black 
hole in flat rather than global AdS space. ^^ 

^*As we will see below, deviations of the black hole from extremality (and deviations of the form of its 
metric from the metric of a flat space black hole) are crucial to dynamics at r — i? ~ i?'', but are small 
perturbations on dynamics when {r — R) ^ R^. 
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In this region it is convenient to work in a rescaled radial coordinate y = -^ and a 
rescaled time coordinate t = -^- Note that the near field region consists of space time 
points with y of order unity (but not too near to unity). Points with y of order -„ (or 
larger) and y — 1 of order 0{E?) (or smaller) are excluded from the considerations of this 
subsection. 

The metric and the gauge field of the background black hole take the form 

dy' 



ds^ =W 



-V{y)dT' + 



V{y) 



+ y'dni 



V{y) 
Ar{y) 



1 



1 






y-l 



+ R^{l + y'^ 



(4.11) 



/^o 



oo 
k=0 



-K ^(0,2fc) 



/^{0,0) =1 
As in the previous subsection we expand 



imid 



iy) = J2^ 



2k irriid 
(1,2A:) 



(y) 



(4.12) 



fc=0 



To determine the unknown functions in this expansion, we must solve the equation 
jj2^rmd _ q^ -^^j^ere D^ is the gauge covariant Laplacian about the background ( 4.11 ). Our 
solutions must match with the far field expansion of the previous subsection, and the near 
field expansion of the next subsection, but are subject to no intrinsic boundary regularity 
requirements. 



At 0{R ) our equations take the form 



1 d f ■^^^ , , d 



where 



y'^ dy\ ° dy 



Vo{y) 



iraid 
{l,2fc) 



{y) 



z>mid 
{1,2k) 



(y) 



(4.13) 



1 

y^ 



and P(i^2k)^y^ ^^ ^ source term determined (recursively) by the perturbative procedure. 
Ignoring the requirements of matching, for a moment, the solution to this equation is 
determined only upto two integration constants at every order. It turns out that 4>7l^2k) (v) 
grows like y (upto possible logarithmic corrections) at large y and grows like , _ .^. as y 
approaches unity. It follows that the expansion ([4.13| ) is good only when 



R^ « (y - 1) « ;^ 



We now work out the explicit solutions at low orders. PJV 
solution for (j)Vl^^\ (y) is particularly simple, and takes the form 

C2 

^(1,0) - ■ 



S(y) 



vanishes, so the 



ci + 



y^ 



1 
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ci and C2 are the two constants. It is easy to check that the matching of 4'7io){y) with 
4''(i*n){^) sets ci = 1. It follows on general grounds that matching with the (as yet unde- 
termined) near field solution forces C2 to vanish. This is because, were C2 to be nonzero, 
it would match onto a near field solution of order O (w) in the near field region (see the 
next subsection for details), violating the requirement that that our solution has a smooth 
R^O hmit. 

We can now iterate the procedure of this subsection to solve to order in R^ in the 
intermediate field region. We find 



imid 
'^(1,2) 



= 1 



^'^'\{y) = -y 



,2-21og(2/2-l) + C3 + ^^ 

y^ -I 



(4.14) 



so that 

Here C3 and C4 are the two integration constants. C3 may immediately be determined 
by matching with the far field solution; it turns out that this procedure also determines 
/i(Q 2) = 0, the chemical potential that was left undetermined in the previous subsection. 
In order to perform this matching we expand <Pi^ (y) about large y 



imid 

h 



(y) = 1 + R' 



-y^ + C3-41og(y) + 0( — 



4„.4^ 



+ OiR'^y 



(4.15) 



The strategy is now to substitute y = -^ in ( 4.15| ) and then to compare with ( [4.10| ). Of 
course one should only compare those terms that are reliable in both expansions. Terms of 



order R^^r"^^ are reliably computed from ( 4.15| ) only when ni + n < 1. Terms of the same 



form are reliably computed from (|4.10| ) only when m < 1. Consequently, the only terms 
that one may reliably compare are those of the form 0{R^r^), 0{R^r'^), 0{R?'r^) together 
with logarithmic corrections. The difference between the sum of the corresponding terms 
(in ( [4.15D and ( |4.1C| )) is given by 



Difi-erence = R^ (^^ + 2^(0,2) log(r) + /x(o,2) - 2 - 41og(i?) 
and vanishes provided /i(o,2) = &nd C3 = —2(1 + 21ogi?) so that 



C3 



imid/ 



r^iy) = l + R'{-y'-2 log(y^ - 1) - 2(1 + 2 log R) + 



C4 



y2_l 



+ 0{R'^) (4.16) 



C4 will be determined below by matching to the near field region. To facilitate this matching 
in the next subsection, we present the expansion of 0™'^(y) expanded around y = 1. 



imid 



R' 



C4 



2(y-i) 



+0{ 



(y-m 



ca 



+ 3 + 2 log 2 + 4 log ii - 2 log(y - 1) + 0(y - 1) 



(4.17) 
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4.2.3 Near Field Region (r - i?) < i? or (y - 1) < 1 

In this subsection we will determine the scalar field in the near field region r — R <^ R. 
More particularly, we will work in terms of a further rescaled radial coordinate z = ^^m-- 
Note that the black hole horizon occurs at z = 0. Note points at finite z are located at 
r — R^R? or y — l^R?. It is also convenient to work with the new time coordinate 
T = Rt = R?T As in the previous subsection, the background gauge field makes a small 
direct contribution to dynamics in this region. However deviation of the black hole metric 
from extremality (and the diff'erence between an AdS and fiat space black hole metric) are 
all important in this region, and have to be dealt with exactly rather than perturbatively. 
In the new coordinates, the metric and gauge field take the form 



l^-<^)-=-w! 



V{z)dT^ + jyy-^ + (1 + R^zfdnl 



■""*"' ='i (' - (TtW ' ° '""''" ^ ""''-'" 



t^o-^R'^''fJ-{o,2k) (4.18) 



k=0 



M(0,0) —1' /^{0,2) — 



i?^ V {R'^z + lf {R^z + lf ^ ' 



= = Az{l + z) + 0{R^) 
As in previous subsections we expand the field (l){z) as 



^T{z) = Y,R'^ct>f,^^,p) (4.19) 



A:=0 

and plug this expansion into the equations of motion. At 0{R'^^) the equations take the 
form 



Ml + z)^y\l,,){z) = Pil,,){z) (4.20) 

where, as usual, P(i2k)('^) ^^ ^ source term whose form is determined from the results of 
perturbation theory at lower orders. We solve the equation ( [4.20D subject to the require- 
ment of regularity at z = 0. It is possible to argue that the solution to 0q 2fe)(-^) behaves 
like z (upto logarithmic corrections) at large z. 

At lowest order (A; = 0) P}iQ\{y) = vanishes, and the unique regular solution for 
(p}iQ)iy) is the constant. Matching determines the value of the constant to be unity. 

At next order, (i.e. 0{R?)) the solution - after imposing the requirement of regularity 
- is given by 

1 (4.21) 

C2)(^) = " - 2 1°S'(^ + 1) - 21og(z + 1) - Li2(-z) 
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where a is the constant which we will now determine by matching with the intermediate 
field solution. Expanding 4>^i{^z) around z = oo we find 



^i^{z) = l + i?2 



7r2 „. ,(\ 



aH 21og2; + C', 

6 \z 



(4.22) 



We now substitute z = ^3- in ( [4. 22]) and then compare with ( [4.17 ). We find a perfect 
match provided a and C4 are chosen to be the following 

C4 = and a = - ( ^ + 3 + 2 log 2 + 8 log i? 

4.3 Perturbation theory at ©(e^) 

We now briefiy outline the procedure used to evaluate the solution at 0{e^). We proceed 
in close imitation to the previous subsection. The main difference is that at this (and all 
even orders) in the e expansion, perturbation theory serves to determine the corrections to 
the functions /, g and A rather than the function (j). The procedure described here applies, 
with minor modifications, to the perturbative construction at 0(e^"*) for all m. 

4.3.1 Far Field Region, r :$> R 

When r 3> i? we expand 



/2°"*W = E^''A°2!2-(^) 



pout 

fc=0 

00 



52°"*(r) = E^''5(2S)W (4-23) 

fe=0 

00 



where 



fc=o 



/-*(r)=y(r), gO-\r) = ^, ^^ W = ^0 ( 1 ^' 



„ -\^p2ifc„ (4.24) 

k=0 
A^(0,0) =1> ^(0,2) = 

As in the previous subsection, we plug this expansion into the equations of motion and 
solve the resultant equations recursively. The equations take the form 

d 



■(r\l + r'fgl^;,,){r) 



I ■ I- . ■ ; .Ti V vfi V ■ / I hource 
dr 



d ( fi2%ir) \ 2(1 + 2r^) „„, , , ^ 
d^ I 1+^-2 ) " r ^^(2 2^0(^) = Source (4.25) 

d ( ,dA"g,,^{r)\ 

-— r ; = hource. 

dr \ dr j 
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and may be thought of as the equations governing sourced linearized fluctuations about 
empty global AdS space with At = I. 

The equations ( [4.25| ) are easily solved by integration. One of the integration constant 
in the first equation is flxed by the requirement that f(2'2k)^^^ ^^ normalizable (see ( |2.11| )). 
The remaining three integration constants (one in the first equation and two in the last) 
will be fixed by matching with the intermediate field solution below. 

The constraints of matching are particularly simple at 0{R^); they require that the 
solutions for g?2o)^ •^(To) ^^^ ^('20) ^^^ ^^^ regular at r = 0. This is because a far field 



solution of the form -c would match onto an intermediate solution of the form —rj^ . But 
this contradicts our basic assumption that our solutions have a smooth i? — t- limit. 
It follows that at 0{R^) ^^ all our functions obey the same equations - and boundary 
conditions - for the 2nd order fluctuations about the supersymmetric soliton and we obtain 
the same (unique) solution 







/(°2!o)(0 



^(2,0) (^) 



1 



At order ©(i?^) we find 



f&ir) 



4r2 (r2 + if 
3 + 5r2 



4(1 +r2) 

1 
'8(1 + r2) 



k 

r2(l + r2)2 



(4.26) 



1 



A out 

(2,2) 



(r) 



l + 2r2 



:i+r2) 



2^2 



+ 



l + r2 
2 



log 



1 + r 



log 



1 + H 



1 + r" 



+ 



+ /ii + 



(4.27) 



Here /c, hi and /i2 are the three undetermined constants, which will be determined by 
matching with the intermediate field solution. To facilitate this determination below we 
end this subsection by presenting an expansion of ( [4.27| ) about r = 



9T\r) 



R' 



+ O(r0) 



fOUt 

J2 



(r) 



-- + 0(r2) 



AT\r) 



1 



+ 0{r^) 



+ i?2 



+ i?^ 



+ 0{R'^ 



+ O(r0) 



4/12 + 1 
4^2 



O(r0) 



+ OiR") 
+ 0{R^) 



(4.28) 



4.3.2 Intermediate field region, r ^ 1 and {r — R) ^ R^ 

As in the previous section, we find it convenient to work with the variables y = -h and 
-5 in the intermediate field region. Recall also that, in these coordinates, the leading 



r 



^^At higher orders the same reasoning does not forbid singularities, but determines them in terms of the 
known intermediate field behaviour at one order lower. 
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order metric has an overall factor of R^. The metric variables that obey simple equations 
have this factor of R^ stripped from them. For that reason we define Here 

r^'iy) = %=9tt and r^^y) = |f = 5rr 

(here g^ji_y are metric components. In a similar fashion we define A^ ~ ~r ~ ^t- With 
these definitions we expand 

oo 

k=0 

oo 

ff2™'(y) = E^''C2%(y) (4-29) 

fe=0 

oo 



k=0 



where 



-'- 1 'JTT) /■} / \ I ^ -^ 



fr'{y)=v{y), 9r'iy) = y^y Ar'{y) = f^o[l-^ 

oo 

„ -X^pSfc,, (4.30) 

fJ'O - 2_^R fJ'{0,2k) ^ ' 

fc=0 
^(0,0) =li ^(0,2) = 

The equations are slightly simpler when rewritten in terms of a new function 

imid / \ 
K^2,2k)iy) = ^0(y)C2%(y) + ^H^ 



where 



Vo{y) = [i-f 



In terms of this function the final set of equations take the form 

dK(^2,2k) jy) _ 
dy 

d ( ^dA^2iy)^y)\ fdK(^2,2k)iy) 



Source 






|;(.V,S,(.))-.%..(.) + ^(^ 



These equations are all easily solved by integration, upto four undetermined integration 
constants (one each from the first and third equation, and two for the second). It will 
turn out that two of these constants are determined by matching with the far field solution 
while the other two are determined by matching with the near field solution. As in the 
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previous section we will find that k order solutions scale like y at large y, but scale like 



2k 



7 — rnr at small z. 



The solution at leading order, R , is given by 



/(^;?)(2/)=«i(i-^ 






aiy 



2ai 2a3 04 

y2 y4 y2 

203^"^ (2a2 - a4)y'^ 



+ 



(y2-l)3 (y2_l)4 (y2_l) 



(4.32) 



as 



Hereai, a^^ 03 and 04 are the four integration constants to be determined by matching. 
Expanding ( 4.32| ) around y = 00 one finds 



fT\y) =«i + ^"^""r""^^ + o f i ) + o(^^^ 



ff2™''(y) 



y vy 

(2q2 — ai — 04) „ / 1 



+ 



+ 0(i?^ 



(4.33) 



^r''(y) =«2 



03 



+ 0(i?2) 



As usual, we substitute as y = -^ and then match relevant terms of (|4.33D and (|4.2g| ) . This 
determines 



02 = — 

/12=- 



«3+4 



(4.34) 



To facilitate matching with the near field region in the next subsection we expand /™'^(y), 5™'^(y) and A'!^^'^{y) 
about y = 1 



ir^y) 



- + 2a3 + Q4 j - (1 + Sas + 2^4) {y-l)+l^ + 20a3 + 3q4 ) (y - 1)^ 



+ o((y-i) 



+ 0(i?2 



gTHy) 



1 + 803 + 4a4 1 + 803 + 804 1 + 803 + 44a4 



+ 



64(y - 1)4 
1 



32(y - 1) 



-n3 



128(y-l) 



y-1 



+ 0(i?2) 



^r'(y) 



+ 03 + 2a3(y - 1) + O ((y - 1)') + 0(fi2 



(4.35) 
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4.3.3 Near Field Region (r - i?) < i? or (y - 1) < 1 

As in the previous section we work with the shifted and rescaled radial coordinate z = ^5- . 
In this coordinate the black hole horizon is at z = 0. As we have seen in the previous section, 
gTT and grr have an overall factor of F? even at leading order. For this reason the natural 
dynamical variables in the problem are 

9TT gtt -I gzz „4 

(here (/^j^ is the metric). For easy of matching with the intermediate field solution however, 
we will continue to use the notation 



fin 

J - gtt 


_ o4 gTT 
-^ X ^2 


5'" = grr 


_ ^ ^9zz 

m i?2 


yl^" = A, 


-R^x^^ 



R 

And so our perturbative expansion takes the form (note the lower limits of the sums) 

00 

k=2 

00 

gT{z)= E^''€2fc)(^) (4.36) 

k=~2 

00 

AT{z) = Y,R"'^il2k)i^) 

k=l 

We now come to an important subtlety of our expansion procedure. First recall that 
the radial coordinate r employed in this paper has geometrical significance; it parametrizes 
the volume of the S^ at that point. For this reason reparametrizations of r do not form 
a symmetry of the equations in this paper, in general. At leading order in the near field 
region, however, the metric metric takes the form 

f = -4.(l + .Mr^ + ^ + .0l (437) 

Note in particular that the size of three sphere (at leading order) is a constant independent 
of z. For this reason the leading order metric equations in the near field region admit 
a whole functions worth (instead of 4 numbers worth) of solutions, parametrized by any 
0{e^R^) redefinition of z coordinate. So without even doing any calculations, we have 
deduced that one linear combination of the three functions is undetermined at leading 
order. 

Now let us move to the next order, 0{I^). As the homogeneous part of the equations 
are same at every order, the same linear combination of second order fiuctuations disappears 
from (i.e. is undetermined by) the second order equations. However the order 'gauge 
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transformation' is now not a symmetry of the 0{E?) equations (because, at this order, we 
see the fact that the size of the sphere is not really constant). So the zero order 'gauge 
transformation function' shows up in the second order equations. As this term comes with 
an explicit B? (without this factor the equations cannot distinguish it from pure gauge) 
it cannot multiply any of the 2nd order unknowns, and so appears as a genuine unknown 
all by itself. The net upshot of all this is that at every order other than the leading, we 
actually do have as many equations as variables. The variables, however, consist of two 
unknown functions at that order coupled with the one unknown 'gauge transformation' at 
the previous order! 

We will now say all of this more precisely. Our equations can be simplified by perform- 
ing the following redefinition of the functions {W is essentially the 'gauge transformation') 

/S2t)(^)=^[1^(l + -il»'(2,2*)(=) 

95..,(.) =c,2,2.,(.) + ^ [-^\ H'<.2.,(.) + y^] i [w,,,,M («8) 

^(2,2fc)(^) =X{2,2k){z) + 2W(2,2k){z) 

In terms of these functions the equations at order B?^ take the form 

— [z{l + z){l + 2zfC{2,2k){z)] = Source 

^ [X{2,2k){z)] - 42(1 + z)C^2,2k){z) = Source (4.39) 

(f 

-^ [W{2,2k-2){z)] = Source 



As we anticipated above, W(2,2k) [z) does not appear in the homogeneous equations at 
0{R^^) as at this order it is pure gauge. But it appears in the homogeneous equations of 
0{R ~^'^). Therefore to completely determine the metric and gauge field (upto integration 
constants) at any given order R^^ , one has to solve one more equation at the order i?2A:+2 
along with all the equations at order R"^^. 

The equations ( [4.39| ) are completely well posed, and may easily be integrated to solve 
for C{2,2k)iz), X{2,2k){z) and W(^2,2k-2)iz) in terms of four integration constants. Two of 
these constants are determined by the requirement that f(2 2k)^^) ^^^ ^?2 2fc~)('^) vanish at 
the horizon z = 0. The remaining two constants are determined by matching with the 
intermediate range solution. 

Solving the first two equations at 0{R^) and the third equation at 0{R^) one can find 
C(2,o)(^)) X{2,o){z) and W{2fl){z) respectively. The solution is the following. 

r ( \-^^^^^_^ 

^(2,0)(^)-^(l + ^)(l + 2z)2 

X(2,o)(^)=^ + A2 (4.40) 
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Regularity at the horizon imphes that 

o 

After imposing the regularity at z = the solution at the leading order 

fi2A)i^) = - ^(1 + 2^) l°g(l + ^) - 3(log2)z + 4/3iz(l + 2z) + 8/322 
€-4)(^) = oo,2nV.^2 [(1 + 2^) l°g(l + ^) + 2^(3 log 2 - 1)] + f'~f' (4.41) 



32z2(l+2)2Lv ; -t>y • ; < V" -b ^J ' 4z(l + z)2 

1 

(2,2) V^; --4 



^f2,2)(^)=-7log(l + ^) + 22/3i 



Here /3i and P2 are the two constants which are to be determined by matching. Expanding 
around z = 00 one finds 

/r(z)=i?4[8/3iz2 + 0(z)]+0(i?6) 

aTiz) =0{R-'^) (4.42) 

Atp{z) =R^ [2Piz + 0{z^)\ + 0{R^) 

After substituting z = ^5- this expansion will match with ( [4.35| ) provided one chooses the 
constants in the following way 

1 
a, = -- 

ai =0 (4.43) 

'^ 8 

In the whole solution at this order there are two constants left undetermined. The first is 
(32 in the near field solution and the second is hi in the far field solution. In the expansion 
of f^{z) the constant (32 appears at 0{R'^z) which is equivalent to a term of O [R'^{y — 1)] 
in expansion of f2^ (y)- Therefore to compute this constant one needs the solution upto 
0{R^) in the intermediate region. Upon determining this solution to 0{R?) one can solve 
for /32 (as well as all the new constants appearing in the 0{R^) intermediate solution) in 
terms of hi. It turns out that 

^2 = l + ^log2-^ 
^ 16 8 ^ 2 

So at the end the full solution to 0{R^) is determined in terms of a single constant, hi, 
which in turn is determined only by the e^ order scalar field analysis. It turns out that 

/ii =0 

5. All Spherically Symmetric Supersymmetric Configurations 

In this section we will analyze the set of spherically symmetric supersymmetric solutions 
of (1.2). The configurations we will find will include the solitons of section ^ (determined 
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more simply than in that section), but will also include several solutions that are singular 
at the origin. In particular (as we have explained in the introduction) we will identify 
a one parameter set of singular supersymmetric solutions which we will conjecture to be 
physical; we will conjecture that these singular configurations may be obtained as the limit 
of nonsingular nonextremal solutions. 

5.1 The Equations of Super symmetry 

The action (|1.2|) is a consistent truncation of AA = 8 gauged supergravity. Hence any 



solution of the equations of motion (2.1C), which saturates the BPS bound m = 3g, cor- 
responds to a supersymmetric solution of A/" = 8 gauged supergravity and consequently of 
IIB SUGRA on AdS^ x S^. Here we present a more direct analysis of the supersymmetry 
equations for the consistent truncation. 
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The supersymmetry conditions in theories of this type have been analyzed in |15, 16 



17| , P?]] . In these works supersymmetric solutions were found for a more general truncation 
of A/" = 8 supergravity to U{1)^ gauged supergravity coupled to 3 hypermultiplets. Our 
theory is a special case of theirs where all three U(l) charges and three hyperscalars are 
taken to be equal. Specializing their results to our theory we find that spherically symmetric 
supersymmetric configurations can be written as follows. The metric, gauge field and scalar 

21 

are 



h? l + p^h? ^ ^ ■^ 



ds = T3, dt^ + ^ ^ .^^^ dp^ + p^h dn^ 

(5.1) 



A = h-^dt, (p = 2v^(/i + p/iV2)2 - 1 

The entire solution is then determined by the single function h{p) which is constrained to 
satisfy the following ordinary differential equation 

{l + p'^h^){3h' + ph") =p\4- {2h + ph'f] h^ (5.2) 

Notice that prime denotes differentiation with respect to the variable p. 

This parametrization of the metric is somewhat different from the one that we used in 
the previous section, so we explain how the two are related. Comparing the coefficient of 
drig in the metric (5.1) to that of (2.9) we see that the two radial coordinates are related 



by 

r' = p'hip) (5.3) 



^"Let us first briefly describe how one could supersymmetrize the action ([1.2|). The bosonic field content 
of our theory is that of minimal gauged supergravity (i.e. the graviton and the (7(1) gauge field) coupled 
to matter (the charged scalar (j>)- The scalar field can be thought of as a member of a hypermultiplet. A 
complete hypermultiplet would contain 2 complex scalar fields. Therefore, to supersymmetrize the action 



(1.2) we have to add one more scalar field, besides the fermions, and the resulting theory is minimal gauged 
supergravity coupled to a hypermultiplet. However, in the set of solutions that we are interested in, the 
additional scalar field can be consistently set to zero so we can ignore it in what follows. 

^^Let us explain our notation in relation to the notation of llq]. We have: rthere ~ Phere, {H\)there ~ 

\H2jthere = [Hsjthere = hhere, Athere = —Ahere, 2 sinh((pthere) = 4>here- 
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Comparing the other coefficients of the two metrics we find 



With these identifications it is a matter of algebra to verify that equation ( ^.71) is sufficient 
for the equations of motion ( 2.1C| ) to be satisfied. 



In summary the most general spherically symmetric supersymmetric solutions to the 
equations of motion ( p.lO| ) is given by the configuration 



fir) 



(2/)/l + p2/j/)2(l+p2/j3) 
1 + />2/l3 



/l2 (5.5) 

^(^) = W) 



0(r) = 2y/{h + ph'/2y-l 
with 



„2 2 



p'Hp) (5.6) 



and h{p) any function that obeys (5.7) 



5.2 Classification of Supersymmetric Solutions 

As we have explained in the previous subsection, supersymmetric solutions to the equations 
of motion are given by solutions to the second order differential equation 

{l + p'^h^){3h' + ph") =p\4- {2h + phy] h^ (5.7) 

In this paper we are only interested in regular normalizable solutions to these equations. It 
is of crucial importance to this section that the condition of normalizability is automatically 
met; an analysis of ( |5.7D at large p immediately reveals that all solutions to this equation 
behave at large p like 

h{p) = l + ^ + ... 

22 ensuring normalizability for all physical fields2^. The importance of this observation is 
the following; one may study the small p behaviour of supersymmetric solutions in a purely 
local manner, without having to worry about when the solutions we study have acceptable 



22t 



Using ([5.3|) and (5.1) we find that this imphes the large r behavior of the gauge field 



so the constant q may be identified with the electric charge of the solution, in the conventions of previous 
sections. 

^^This fact has a natural explanation from the viewpoint of the dual A/" = 4 Yang Mills field theory; 
a deformation of the Lagrangian of that theory by only mass term TrX^ + TrY^ + TrZ^ preserves no 
supersymmetry. 
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large p behaviour, as that is always guaranteed. This fact allows us, in this section, to 
use local analysis to present a simple classification of normalizable supersynimetric solu- 
tions. Relatedly, supersymmetric solutions may be obtained by solving ( |5.7D as an initial 
value problem with initial conditions set at small p. This is numerically and conceptually 
simpler than the boundary value problem we would have to solve by shooting methods off 
supersymmetry. 

It remains to impose the condition of 'regularity'. Let us first explain what we mean 
by this term. We call a supersymmetric configuration 'regular' if it can be regarded as 
the limit of a one parameter set of smooth nonsupersymmetric (and so non extremal) 
solutions to the equations of motion ( |2.10| ). While every smooth supersymmetric solution 
is automatically 'regular', a singular susy solution may also be 'regular', if its singularity 
can be removed upon heating the solution up infinitesimally. 

Solutions to ( ^??| ) can develop singularities only at /? = 0. In this subsection we will 
classify all possible behaviours of solutions to ( ^.7D near p = 0. In a later subsection we 
will then go on to present conjectures about which of these solutions are 'regular'. 

In order to investigate possible behaviours of solutions to ( ^.7[ ) at small p we plug in 
the ansatz h{p) = —^ into the equation. It is easy to check that the only values of a that 
solve the equation near p = Q are a = 0, |,1,2. It is also possible to demonstrate (see 
below) that there is a unique solution with a = g . On the other hand solutions with a = 
and a = 1 both appear in a one parameter family. Finally, the generic solution to the 
differential equation has a = 2; solutions with a = 2 appear in a 2 parameter family. ^^ 

5.2.1 h{p) ^ p-i 

As we have mentioned above, there is a unique solution with a = g. This solution may be 
expanded at small p as follows 

h{p) = p-^l^ + - p2/3 _ JIL p2 ^ 0(8/3) (5 8) 

^^' ^ 2Q^ 20956'^ ^^ ' ^ ' 



We now present a crude estimate for validity domain of the expansion ( ^.^ ) . Note that the 
formal procedure that generates the series expansion ( |5.8| ) treats the term proportional to 
4 (in the RHS of ( ^.7D ) as subleading to the term proportional to h. This procedure is 
valid whenever p ^ 1; as a consequence we expect the expansion ( |5.8D to be valid whenever 
p <^\ but to break down at larger values of p. 

For /9 ^ 1 the metric, gauge field and scalar corresponding to this solution take the 
following form 

8 ^ 

A{r) ^ rdt (5.9) 

4 
6r 



^^One special exact solution of ( p.7| ) is h(p) = 1 + -%, the so-called "superstar" [p8|. This solution has 
a = 2. For this solution we notice that the scalar field is not turned on (0 — 0). So in a sense it is 
qualitatively difi'erent from the "hairy" configurations of interest to us in this paper. 
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where we used the relations (|5.5D and (|5.q ) to bring the solution in the form of (^ 

We will denote the distinguished singular solution of this subsection by S. We will 
now explain that there is a sense in which 5" is a fixed point of the equation (|5.7|) viewed as 



a flow equation in the variable logp = x (see |29, 3C, 31 1 for similar discussions in distinct 



but similar contexts). For this purpose we redefine 

Kp) = e-hfix) 

The differential equation ( |5.7D becomes 

9r(l + f) + 3/'(2 + 3/2/' + Wf) + 8(/^ - /) - 36e^"/V' = 

For very small p (that is for x — )• — oo) the last term in the equation can be ignored, and the 
equation becomes approximately time translation invariant (or an autonomous equation, 
in the language of dynamical systems). With this approximation the system has an exact 
solution f{x) = 1, which is precisely the leading small p approximation to the solution S. 
We will restrict attention to large negative values of x in the rest of this subsubsection, 
and so study the truncated equation 

9/"(l + f) + 3/(2 + 3/2/' + 10/=^) + 8(/4 - /) = (5.10) 

Let us consider a small perturbation about / = 1, i.e. we set 

fix) = l + egix). (5.11) 

To linear order in e the ( ^.10 ) turns into the linear ODE 

3g" + 6g' + Ag = (5.12) 

The two linearly independent solutions to this equation are 

9{x) = e"^ 



(5.13) 
Note that the real part of the each of these eigenvalues is negative, which demonstrates 



\/3 



that / = 1 is a stable fixed point of the dynamical system (5.10). Moving back to the 



variable p, it follows that arbitrary small perturbations around the solution S behaves like 

h(p) Rs p"2/3 M ^ (J _ cos I ^ log /9 + a 

V P Vv3 

Note that all perturbations die out for p ^ e (this is a restatement of the fact that / = 1 
is a stable fixed point). 
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5.2.2 h{p)^ho + 0{p'^ 



We now turn to regular solutions to (5.7), i.e. solutions with a = 0. Such solutions appear 



in a one parameter set, labeled by /iq = h{0). The small p expansion of ( |5.7D is given by 
Kp) = ho + l{hl- hDp" + \{hl- bhl + ^hl)p^ + 0{p') (5.14) 

The solutions of this subsection are simply the solitons studied in section |3[ These solutions 
were generated perturbatively (i.e. at /iq — 1 small) in section ^. 

Let us now turn to the opposite limit of large /iq. Let us first inquire as to the 
validity domain of the expansion ( 5.14 ). As the term h^p'^ on the LHS of (5/7) is of 



order p"^, the formal process that generates the power series expansion ( ^.14[ ) treats this 
term as subleading compared to unity. This is actually correct only when h^p'^ ~ 1. It 
follows that when /iq is large the series expansion ( ^.14 ) will break down at the small value 

_3 

p ~ Prb ~ ho ^ 



_ 3 

When ho is large, the expansion ( p. 14 ) does not apply in the range p ^ ho^ ■ If p ^ 1 



in this range, however, the general arguments presented above guarantee that our solution 
behaves like p~°' for one of the allowed values of a described above. What solution does the 



expansion ( 5.14 ) match onto in this range? A clue to the answer to this question is given by 
noting that the value of h{p), at the point of break down of ( |5.14| ) is approximately given 

_ 2 

by ho ~ p^i^^ . Thus the solution ( 5.14| ) could smoothly match onto the special solution S of 



previous subsubsection, at p ~ prb- This suggests that the special solution S of the previous 
subsubsection is the limit as /iq — )• oo of the solutions of this subsubsection. We now present 
numerical evidence that strongly supports this guess. In Fig. (^ we present numerically 
generated plots of the regular solution parametrized by ho for successively increasing values 
of ^0- Note that, as ho increases, the entire profile of the solution approaches a limiting 
shape, with a sharp spike near p = 0. The spike becomes sharper as we increase /iq, while 
the solution at larger values of p remains almost unchanged. The limiting solutions indeed 
appears to be the special solution S of the previous subsubsection (denoted by the solid 
line in Fig. |3|). Thus it appears that the solution S forms the end point of the family of 
regular super symmetric solitons. 

We will now study in more detail how solutions with large /iq approach the special 



solution S. We work in the language of the dynamical system (5.1C). We are given a 
function /(x) that starts out, at large negative values of x (small p) as 

f{x) = hoe"-^{l + 0{e^n). (5.15) 

We wish to study how f{x) evolves under (|5.10D at later times. We are interested in the 
limit in which ho is large; as we have argued above, we expect f{x) to increase to a value 
of order unity at a time xq ~ — |ln/io, and thereafter stabilize exponentially to the fixed 
point / = 1. Note that xq <^ —1 (this follows because of our assumption that ho is large) , 
so that / should settle down to very near unity well within the domain of applicability of 
the dynamical system ( ^.1C1| ). ■^^ 

^^In this language, the conjecture of the previous paragraph is equivalent to the assumption that this 
fluctuation lies within the domain of attraction of the fixed point f — 1. 



35 



h{p) 




— q — ^ 



P 

Figure 3: Convergence of the numerical solutions 
for the regular solitons to the special singular so- 
lution S as we increase Hq = h{Q). The black line 
corresponds to the solution S with p""^!^ behavior 
near p = 0. The blue lines correspond to regular 
solitons of Kq = 2, 3, 4, 6, 8, starting from the low- 
est blue curve and going up. 




log/io 
Figure 4: The damped oscillations of q around 
the critical value qc for large /iq. 



Let us now compute /(xi) for some fixed (/iq independent) x\ that obeys^^ 

xo ^ xi ^ —1. 



In order to do this we recall that the equation ( |5.10| ) is invariant under translations in x. 
Now the judiciously chosen translation 



x^ 



31n/i( 







eliminates the /iq dependence of the initial condition ( |5.15| ) . Let x^"^) be the (unique, /iq 

2i' 



independent) solution to ( 5.1C| ) that reduces at early (large negative) times to x(x') = e s . 
It follows that the solution of interest to us is 



31n/i| 



On 



/(x) = x(a: 

The key assumption of this section, is that the function x li^s within the domain of 
attraction of the fixed point / = 1 (as we have seen above there is impressive numerical 
evidence for this assumption). If this is the case it follows from ( 5.11 ) and ( ^.13| ) that at 
large x' (i.e. for x 3> xq;) 

x' 
^[x') = 1 -h ^e"^ cos(^ -h S) 

for some unknown, order unity constants A and b. It follows that 



V3 



/(xi) ^ 1 + Ae-2 '"S'^o+^i cos ^ log(/io + xi) + (5 (5.16) 

Of course the xi dependence of this result may be absorbed into a redefinition of A and b. 



Recall that (5.10) is valid only for large and negative x\. 
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Let pi = e^^ . ( 5.16| ) gives us a formula for h{pi) and h'{pi) for the solution of interest; 



these values may be used as an 'initial conditions' to generate h{p) for all p > pi . The 
resultant solution will take the form 

h{p) « hs{p) + 5h{p) 



where hs{p) is the special solution S and 5h{p) is a small fluctuation (of order ~ C'(— 3-) 

) about this solution. To leading order in this small parameter, the function 5h{p) obeys 
a linear differential equation, and so depends linearly on /i(pi) and h'{pi). It follows that 



(5.16) then determines the behaviour of every observable (like the charge) of the solution 
that depends only on the behaviour of h{p) for p of order unity or greater. In particular it 
follows that the dependence of the charge of solutions on ho is given approximately by 



3 !„„;,„ / \/3 . 

Y 



g(/io)~9c + ^e-t^°s/iO(,os [l^logho + 6] (5.17) 



for some constants A, 6. While A and 6 can only be determined numerically, we have a 
sharp analytic prediction for the form ( ^.17 ). A similar formula applies for the vacuum 



expectation value of the operator dual to i;^ as a function of /iq- 

We have verified the prediction ( |5.17 ) numerically; in Fig. ^ we present a plot of the 



rescaled oscillations of q about Qc- This graph displays precisely the damping (reflected in 
the ho dependent renormalization of the y axis in Fig. ^) and the oscillations predicted by 
(5.17). We will give more details below. 



5.2.3 h{p) « ^ 

Next we consider solutions with a = 1, i.e solutions that behave near p = like ^. The 
one parameter set of these solutions may be labeled by o. At small p our solution takes the 
form ^P{p) where P{p) is a regular power series. The first few terms in the power series 
expansion are given by 

w s ^ f 1 18a4-5 2 -90a^ + 31 o ,^, .A 



The formal procedure that generates the power series (5.18) treats the term proportional 
to unity (on the LHS of (5.7)) as subleading compared to the term proportional to h^p^. 
When a <C 1, this is justified only when p <C a?- Consequently we expect the expansion 



(5.18) to break down at pg^ ~ a?. 



As in the previous subsubsection, the solution presented here must reduce to one of the 

other solutions of this section when a? -^ p -^ 1. Noting that, at the point of breakdown 

_ 2 

of (5.18), the function h may be estimated by /i ~ ^ ~ p ,^ , it is natural to guess that the 



solution of this subsubsection tends to the special solution S in the limit of small a. We 
now present strong numerical evidence in support of this guess. In figure (g) we present 
numerical plots of the solution of this subsection for a range of decreasing values of a. Note 
that the solution converges to the solution S (denoted by the solid line in Fig. (|5|)) at small 
a. 
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Figure 5: Convergence of the numerical solutions 
for singular solitons with an - singularity to the 
special singular solution S as we decrease a. The 
black line corresponds to the solution S with p"^/^ 
behavior near p = Q. The red lines correspond to 
singular solitons of a = 0.35, 0.5, 1, 2.5, 5, starting 
from the lowest red curve and going up. 
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Figure 6: The damped oscillations of q around 
the critical value q for small a. 



For p <C 1 the metric, gauge field and scalar corresponding to this solution take the 
following form 



ds^ w -r'^dt^ + -^dr"^ 
a* 

A{r) « ^dt 



r'^d^l 



(5.19) 



Precisely as in the previous subsubsection, we can analytically study the approach of 
the solution with small a to the solution S. Repeating an analysis very similar to that of 
the previous subsubsection, we conclude that the dependence of, for instance, the charge 
of the solution on a is given by the formula 



g(/io) ^qc + Ae^^°^°- cos {V2,\og a + 5\ 



for some constants A, 5 which cannot be determined analytically. We have verified this 
prediction numerically (see Fig. ^. A similar formula applies for the vacuum expectation 
value of the operator dual to the field (p^ as a function of /iq . 

5.2.4 The generic solution, a = 2 

Finally we move to the case m = 2. Now we find a two parameter set of solutions, labeled 
by two arbitrary constants a, 6, which have the form 



h{p) 



1 



a + 



1 



2a 



■P 



(62 - l)(3a(562 



1) 



24a4 



/ + 0(p 



12\ 



(5.20) 



The value 6 = 1 is special; as we have already remarked h = 1 + -% is an exact solution; 



hints of this fact are already visible in the expansion ( ^.20 ). It follows in particular that 
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the values 6 = 1 lies outside the basin of attraction of the fixed point S at least when b = 1. 
A very rudimentary numerical investigation suggests that this is also true for all values 
of a at (for example) 6 = 2. Although we have not carefully investigated this question, 
it seems possible that the solutions with a = 2 are completely disconnected from all the 
other solutions studied above. 

For the reasons outlined in the previous paragraph, the 'generic' solutions of this 
subsubsection will make no further appearance in our paper. We suspect that all - or at 
least most - of these solutions are genuinely singular, in the sense that they cannot be 
regarded as the limit of smooth solutions. We leave a fuller study of these solutions, and 
their possible physical significance, to future work. 

5.3 'Regular' supersymmetric Solutions 

In this section we present the results of a numerical analysis of the space of 'regular' 
supersymmetric solutions. Let us first describe what we believe the space of these solutions 
to be. The smooth supersymmetric solitons (with a = 0) of the previous subsection are 
clearly regular. However this space of solutions ends at finite charge (see below) as ho — )• oo. 
As we have described above, this line of solutions spirals into (and ends in ) the special 
solution S. We have also seen above that another line of solutions - those with a = 1 
and small values of a - spiral out of the solution S. As we will see below, the charge 
of this new line of solutions increases without bound at large a. It thus seems that the 
solutions with a = and a = 1 may be regarded as two different segments of a single 
line of supersymmetric solutions. The two segments are joint together (by a very intricate 
non intersecting double spiral structure) at the special solution S. At least one member of 
this line of solutions exists at every value of the charge, and so constitutes a candidate end 
point of the phase diagram Fig. g. We conjecture that it is indeed the case that hairy black 
holes at every value of the charge exist for all energies above the BPS bound. In the BPS 
limit, these solutions reduce to some configuration on this special line of solutions; either 
to the smooth supersymmetric soliton (at small charges) or the a = 1 solutions (at large 
charges). In particular we conjecture that all solutions on the special line described in this 
are 'regular', where this word is used in the sense specified in the previous subsection. We 
now proceed to study these conjecturally regular solutions in more detail. 

5.3.1 Solitons 

We first present the numerical analysis of regular solutions of ( ^.71 ). For this we fix the 
initial conditions /i(0) = /iq ^ 1 ^^ and h'{0) = and integrate the equation outwards. 
For each value of /iq we compute the solution h(p) numerically and we evaluate the charge 
(/(/iq). The results can be seen in Figs 0, ^ and |9|, in various magnifications. ^^ 



''The condition /lo > 1 is necessary since tlie scalar field at p = is given by 0(0) — '2\/hQ — 1 and has 



to be real by assumption (2.9). 

^*To partly check the validity of the numerics it is easy to perform a perturbative analysis of equation 
( |5.7| ), similar to that of the previous sections i.e. in a small amplitude of h{p) — 1. One finds agreement 
between numerics and perturbation theory (i.e. convergence of the perturbative solution to the numerical 
one, for small enough values of ho — 1). Notice that the case /lo = 1 is precisely empty AdS. 
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Figure 7: Charge q of spherically symmetric su- 
persymmetric regular solitons as a ftmction of the 
value ho = h{0). 
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Figure 8: The same graph with different scales 
on the axes, where we can see the maximum 
charge. 




log/io 
Figure 9: Magnification of the previous graph. We wee the next oscillation around Qc. 

The most striking feature of the numerical analysis is the existence of a maximum 
value^^ 

Qmax ~ 0.2643 (5.21) 

for the charge of regular supersymmetric solitons. This charge is obtained for the value 
^qmax ^ 9.821 of the initial conditions at the center. The existence of a maximum charge 
for these solitons was also noticed in |1^]^'^. For higher values of /iq the charge of the 
solution starts to decrease and asymptotically it approaches the limiting value 



0.2613 



(5.22) 



as ho — )• oo. 

A more careful analysis reveals that the convergence of the function q{ho) towards the 
critical value Qc is not monotonic. Instead, the function q{ho) undergoes slow oscillations 
around the critical value, as can be seen in Fig. |^ and in more detail in Fig. |^. These 
oscillations are periodic, with damped amplitude, if expressed in terms of x = log/io- As 



we explained in section 5.2.2 the asymptotic form of these oscillations can be determined 



analytically by matching the regular solution for large /iq to the special p ^' ^ solution S 



^^We have solved the equations numerically using Mathematica. 

''''Notice that we are working in sliglitly different conventions from ||lq], in which qthere = 2qhere- This is 
consistent with the maximum charge {qm)there = 0.529 reported in that paper. 
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and we have the following asymptotic formula for large /iq 

q{ho) ^ Qc + Ae~^^°^''" cos{ujlogho + 5) (5.23) 



In section 5.2.2| we saw that the period of the oscillations and the damping constant can 



be determined analytically from the matching procedure to be 

7 = 3/2, oj = V3/2. (5.24) 

while A and 6 cannot be fixed in this way. If one tries to fit this formula to the numerical 
data one finds 

7 w 1.50, WW 0.87, A^-0A9, 5 ss 0.12 (5.25) 

which are in very good agreement with the exact values. 

The number of regular solitonic solutions as a function of the charge are as follows: for 
small enough charge there is only one solution. As we increase the charge, at some point 
we hit the first oscillation around Qc, which increases the number of solutions to three. 
Increasing q further we encounter the second oscillation and we have five solutions, and 
so on. As we approach the critical value Qc from below the number of solutions is always 
an odd integer which goes to infinity. Now let us consider the large charge behavior. For 
Q > Qmax we have no regular solitonic solution. As we decrease the charge and we go below 
Qmax we first find two solution. As we decrease further we encounter the first oscillation 
above qc, giving us four solutions, then the second oscillation to six solution and so on. 
Hence for q > qc we always have an even number of solutions (possibly zero) which tends 
to infinity as we approach qc from above. 

Notice that since the BPS bound m = 3q is satisfied for all of these solutions, the 
figures |^^ also show the dependence of the mass of the solution on the value of the field 
at the center. This qualitative behavior, i.e. the existence of a maximum mass, and of a 
slightly lower critical value of the mass which is approached asymptotically for large central 
density via a function which undergoes damped oscillations, is typical in related problems 
in general relativity [ |32| , 33, 34]. ^^ To our knowledge, however, this is the first time such 
behaviour has been observed in family of supersymmetric solutions. 

Let us now study the expectation value of the scalar operator dual to (j), which we 
denote by (O,^), as a function of ho. Since C^ is an operator of dimension A = 2, its 
expectation value can be determined from the large r expansion of (p as 



^^ Generally, when a family of gravitational solutions has the property that their mass has a local maximum 
for some value of the initial conditions at the center, it is the sign that one of the two branches (to the left 
or right of the local maximum) is unstable (under radial perturbations) and thus unphysical. This is the 
analogue of the "Chandrasekhar instability" : if we expand the equations of motion around the solution at 
the local maximum of the mass they have a zero mode, since the total mass does not change to first order 
in the perturbation. Generically this zero mode will be stable on one side and unstable (i.e. tachyonic) on 
the other side of the local maximum. This is what happens for example in the case of boson stars |p5|] . In 
our case the solutions are supersymmetric for all values of ho- It would be very interesting to check what 
this implies about their stability. 
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We plot the results in figures 16,17,1^ in appendix |C|. The qualitative behavior is similar 
to that of the charge q: the expectation value (O^) is an increasing function of /iq up to a 
maximum value 

{0^)max ~ 1.8906 (5.26) 

which is realized for ho ~ 6.580 and then decreases and approaches the asymptotic value 

{0^)c w 1.8710 (5.27) 

as /iQ — )• oo, while performing small oscillations around it. Again the oscillations can be 



determined following the logic of section 5.2.2 and are captured by the formula 



{O^){ho) ~ {0^)c + Ae-^'''^'''cos{ojlogho + 5) (5.28) 

The analytic prediction is 7 = | , uj = ^. If one tries to fit this formula to the numerical 
data one finds 

7 Pi 1.50, w«0.87, A ^-0.66, 6^0A8 (5.29) 

which are in good agreement with the exact values. 

Before we proceed let us point out that the value of Hq at which we have the maximum 
charge {hq^^^ ~ 9.821) differs from the one at which we have the largest expectation value 
{Off)) max-, which turns out to be h(p a^^^ ~ 6.580. More generally, and as we will see more 



clearly in subsection 5^ , while there are pairs of regular solitonic solutions with the same 
charge q or the same expectation value {O^), there are no such pairs which have the same 
q and same (O^) simultaneously. 

5.3.2 Solutions with a = 1 

We now present the results of a numerical investigation of the other segment of the line of 
(conjecturally) 'regular' supersymmetric solutions: those whose small p behaviour is given 
by a/p for small p. We compute the entire solution numerically and calculate the charge 



g as a function of a. The results are shown in Figs 10,11|,1^. As we see in the figures 
the charge of this family starts (at small a) precisely at the point q = qc ( ^-221) where the 
family of regular solitons ended, then as we increase a the charge seems to decreases, down 
to a minimum value 

qmin = 0.2605 (5.30) 

and then increases all the way to arbitrarily large values. 

As before, a more careful analysis shows that the entire radial profile of solutions with 
- singularities converges to the special solution S in the limit a — ?■ 0, as shown in Fig. |5|. 
Again, a closer inspection shows that in the regime between a = and the point where 
Q = Qmin the function q{a) is not monotonically decreasing, but rather is undergoing small 
damped oscillations around the value qc as a function of log a. This is shown in Fig. |6|. 
For small values of a the form of these oscillations can be determined by the matching 
procedure discussed in section ^.2.3 and we find the following formula 



qia) ^qc + ^e^'°S" cos(w log a + 5) (5.31) 
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Figure 10: Charge q of spherically symmetric 
supersymmetric solitons with a singularity of the 
form - at p = 0. 
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Figure 11: Detail of the previous graph with 
different scales on the axes, where we can see the 
minimum of the charge near a = 0. 
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Figure 12: Magnification of the previous graph. We wee the next oscillation around Qc 



with the analytically determined values (see 5.2.3| ) 7 = 3, uj = ^/3. From the numerics we 
find 

7 ^ 3.00, oj pa 1.73, A ^ 0.18, 6 w 0.70 (5.32) 

which are in good agreement with the exact values. 



We find similar behavior for the expectation value (O^) as shown in Figs 20,21 



m 



appendix ^: the expectation value starts at the point ( |5.27D where the regular soliton 
family ended, it then goes down to 



{0^)r 



1.8658 



(5.33) 



and then increases indefinitely. 

Finally we have the following oscillatory behavior for small a which is shown in Fig. 

{0^)ia) « (0^), + ^e^'°s"cos(wloga + 5) (5.34) 

with the exact values 7 = 3, uj = v3- From the numerics we find 



7 ^ 2.97, oj « 1.74, A ^ 0.59, 



-0.34 



(5.35) 



Let us mention that the numerical results depicted in Fig. 14 agree with the perturba- 
tive analysis of section |^ in the regime of small q. According to the results of that section 
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we expect that for small q the expectation value {O^f,) goes like 

(0^)=4^+... 



One can indeed verify that the small q behavior of the curve in Fig. 14 agrees with this 
result. 



5.3.3 An analytic solution at large charge 

While we have no analytic control of a = 1 solutions in general, we can see from the 
numerical analysis that for large a the charge of the solution can be well approximated by 

2 

the formula q = ^ + subleading. In fact, in the limit of large a one can find an analytic 
form of the solution as follows: let us consider the first factor (1 + p^h?) on the LHS of 
equation ( |5.7| ). At small value of p the term p^hp' dominates over the 1 since by assumption 

Hence it is not unreasonable 



/i ~ - . At large values of p the same is true since h 



1 + 7?- 



2<? 
7 



to assume that in the limit of very large a we can make the approximation 1 + p^/i^ ~ p^h^ 
for the entire range of p. Then the differential equation (| 



becomes 



phi^h' + ph")= A-(2h + ph' 



(5.36) 



This equation can be solved exactly^^ and if we impose the desired behavior near p 
the solution is 



hooip) 



1 + 



(5.37) 



It is not hard to check that in the large charge limit the numerical solutions do indeed 
converge towards the solution ( ^.37 ) in the entire range of p, as shown in figure 13. The 



solution (|5.37 ) goes like - near p = and like 1 + j-j for large p. As we said this implies 



that the charge q ^ \- One also finds that in this limit the expectation value goes like 
{0(j)) ~ a^. So in the limit of large q we have 

{0^)=Aq + ... 



which describes the behavior of the red curve in Fig. 14 for large q 



5.4 Phase Structure of 'regular' supersymmetric solutions 

Let us now put everything together and describe the space of supersymmetric solutions. 
In Figs 14,1^ we show the expectation value {O^j,) and the charge q of the family of regular 
solitons (blue curve) and those with an ^ singularity (red curve). As we explained above 
the two families meet at the solution S with p~'^'^ behavior, which is denoted by a black 
dot. Near the point S the two curves develop into two intertwined spirals which are 
asymptotically described by equations ( ^.23 ), ( 5. 28 ), ( 5. 31 ), ( 5. 34 ). We zoom into the point 
S in Figs. UJ2^ in appendix 0. 



^The general solution of (|5.3(^) is h{p) — . 1 
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Figure 13: Convergence of the numerical solutions for singular solitons with an - singularity to 
the family of approximate solutions /loo(p) ^ we increase a. We plot the ratio of the two functions 
for various values of a and we see that it converges to 1 as we raise a. The values plotted are 
a = 0.1, 0.2, 0.5, 1.5, from top to bottom. 
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Figure 14: Expectation value (O^) vs charge q 
for the family of regular solitons (blue) and for the 
family of solitons with an - singularity (red) . The 
two curves meet at the point denoted by the black 
dot which corresponds to the special solution S 
with p"'^!'^ behavior. 
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Figure 15: Detail of the previous graph around 
the point S where the two families meet. The 
blue curve is the regular soliton and the red curve 
is the soliton with the - singuality. 



From these figures we see that the curves are non intersecting, which means that if 
we fix the charge q and the expectation value (O^) there is at most one solution. If we 
consider the number of solutions as a function of the charge q leaving (C?^) arbitrary we 
have the following pattern: for small enough q we have only one regular solitonic solution. 
As we increase q we first encounter the point qmin where two new solutions appear, bringing 
the total number of solutions to 3. Increasing q we cross a point where two more regular 
solutions are added and the total number of solutions becomes 5. This pattern continues as 
we approach the point qc with an ever-increasing number of total solutions (at each step we 
add, alternatively, either two regular or two singular solutions). Notice that this number is 
always odd. After we cross the point qc the pattern is reversed, we successively lose pairs of 
solutions until we end up with a single singular solution for q > qmax- Similar statements 
hold if we look for solutions with given expectation (C^) instead of given charge. However, 
the fact that the curves in these figures are non-intersecting, means that if we specify both 
q and (O^ there is always at most one solitonic solution with these values. 
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The phase diagram we have proposed for our gravitational system is depicted in Fig. 
|2|. In this diagram we have included a phase transition curve that meets the BPS line near 
q = q^K, 0.2613; we will now explain our rational for doing so. In this paragraph we assume 
that both the regular and the singular super symmetric solutions found before, may each be 
obtained as a limit of non singular non extremal hairy black hole solutions. As we saw above 
for any given value of q we may have either one, or a larger odd number of supersymmetric 
configurations depending on whether q lies outside or inside the interval {qmin,qmax) = 
(0.2605, 2643). It follows that there may exist more than one near supersymmetric regular 
hairy black hole solutions in the charge range q S {qmimQmax) approximately centered 
around qc, on which we now focus. These configurations differ by the expectation value 
(O^) of the operator dual to the field (p. Let Sji{q, 6e) denote the entropy of the hairy black 
hole that reduces to the regular soliton with the largest value of (O^) when 5e — )• (here 
5e denotes the energy above BPS). Let Ssiq, Se) denote the entropy of the hairy black hole 
that reduces to the singular supersymmetric solutions with the smallest value of (O^). We 
suspect that 

SR{q,6e)> Ss{q,6e) when q < qp{Se) 
Ss{q,Se) > SR{q,6e) when q > qp{6e) 

for some qp{Se) such that liuiSe~^o qp{6e) £ [qmimQmax)- Moreover we suspect that the 
entropies of the plethora of intermediate phases that appear at charges near to qc are 
always smaller than either Ssiq,6e) or SR{q,Se). In other words we suspect that our 
system undergoes the micro canonical analogue of a a single first order phase transition at 
q = qp{6e); this is the black curve we have depicted in Fig. |2[ The phase transition curve, 
which originates at the BPS line, could either extend all the way to the phase transition 
curve between RNAdS and hairy black holes, or could terminate somewhere in the bulk 
of the hairy black hole phase, at a triple point analogous to the water steam system. Of 
course the considerations of this paragraph have been highly speculative. It would be very 
interesting to investigate this further. 

Before we continue we would like to mention that it would be important to clarify the 
stability of these solutions under linearized perturbations. As mentioned in footnote ^ 
on general grounds one would expect regular solutions past qmax to be unstable. On the 
other hand our solutions are supersymmetric and from this point of view it would seem 
more natural to believe that they are stable. This is a confusing issue that deserves further 
study. 

6. Thermodynamics in the Micro Canonical Ensemble 

In this section we present thermodynamical formulae for RNAdS black holes, the super- 
symmetric solitons, and hairy black holes, in a small charge and near extremal limit. We 
also demonstrate that the leading order thermodynamical formulae for hairy black holes 
are reproduced by modeling them by a non interacting mix of a soliton and an RNAdS 
black hole with ^u = 1. 
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6.1 RNAdS Black Hole 

The basic thermodynamics for an RNAdS black hole is summarised by the following for- 
mulae 

m 



(6.1) 



M 

- Ar2 


= ^i?2[l + i?2 + ^2] 
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- Ar2 


-1^' 


S 




- Ar2 


= ttR^ 


T 


-..^[1 + 2^' ^' 



were Q is the charge, M is the mass of the black hole, S is its entropy, T its temperature 
and /U its chemical potential. Note that 

^2<(i + 2i?2). (6.2) 

(this follows from the requirement that R is the outer rather than the inner event horizon 
of the black hole) . 

In this paper we are interested in small RNAdS black holes - i.e. black holes with 
m <C 1 and q <^ 1 that are also very near extremality. The mass of RNAdS black holes 
at fixed charge is bounded from below by the mass of the extremal black hole of the same 
charge; at small q we have 

m>3{q + q'^- 2q^) + ©(g^) (6.3) 

For every pair {m, q) that obeys this inequality, there exists a unique black hole solution. 
In this paper we are interested in black holes whose mass above extremality of of order 
0{q^). ^^ For this reason we define the shifted and rescaled mass variable 5"^ by 

,52^3 = m-2,{q + q^ - 2q^) (6.4) 

We are interested in g ^ 1 but 6 of order unity. In this regime the entropy, chemical 
potential and temperature of this black hole is given by 



3 
7rg2 



2\/2 + hV^S - qV2\ q+{^=- 10^3(5 + 33^2 - ^-y^ \ q^ + O {q 




,= 1+ 2-JisU+r-:r-6+^-^ U' + Oi,') 



ttT 



6 (3V25^ - 10V36^ + 24^/3) q , ,, 

Tz M ^+°''' 



(6.5) 



^^Note that the mass of an extremal black hole, at charge g, exceeds the mass of a BPS black hole at the 
same charge by 3q^ — 6g^ + 0(g''). Consequently the deviation of the mass of our black holes from the BPS 
bound is given by 3g^ + ((5 — 6)g'', and in particular is 0{c^) rather than ©(g'"*). 
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Although the black holes we study are very small, their temperature is very small (it 
scales like y^) because we focus on the near extremal limit ^^. Moreover the black hole 
temperature decreases as we decrease 6, reaching zero at 6 = + 0{q). 

Note also that the chemical potential /x of these black holes is unity at leading order. 
The first correction to this leading order result is of order 0{q) and is positive when ^■^ < 6 
but negative otherwise. This already suggests that RN-AdS black holes with 6^ < 6 are 
unstable to super radiant decay; we will see below that this is indeed the case. As we will 
see below, the end point of the resultant tachyon condensation process is a hairy black 
hole. 

Finally note that the radius of the black holes we study is O(y^) so that the entropy 

3 

is of order 0{q2). 

6.2 Supersymmetric Soliton 

The mass and charge of the supersymmetric soliton are given by 



3 /e2 £4 e6 igggS 

m =- I ! I I h O (e 

4 V 4 192 1920 2211840 ^ ' 

^ 2 \ 8 384 3840 4423680 ^ ' 



(6.6) 



Note that 



(6.7) 



m = 3q 
Of course the soliton is dynamically stable as it is supersymmetric. It carries no entropy. 



6.3 A non interacting mix of the black hole and soliton 

In this subsection we will determine the thermodynamics of a hypothetical non interacting 
mixture of the small black holes and the supersymmetric solitons of the previous subsection. 
Of the net mass m and charge q of the system, let mass 3qs and and charge qs lie in 
condensate so that the mass and charge of the black hole are given by 

rrih = m — 3qs 

(6.8) 

qb = q-qs 

The charge qs is determined by maximising the entropy of the system, which determines 
the black hole chemical potential to be unity. This condition gives 

mb = 3qb + 3ql (6.9) 

^*The reader may worry that the appearance of inverse powers of 5 in ( ^.5| ) signify that black hole 
thermodynamics degenerate in the extremal limit; this, however, is not the case. At extremality 5 = 5ext{q)- 
The function 5sxt{q) starts out at 0{q) and so is small at small charge, but does not identically vanish. 



Infact Snxtig) may be determined as a power expansion in q by equating the temperature in (6.5) to zero. 
Plugging this function into the remaining expressions in (pj) , yields nonsingular, analytic expressions as a 
function of q. 

^^In contrast small black holes in [[L4| all had a very high temperatures. 



48 



(this exact formula may also be verified to 0{q'^) by setting 6^ = 6 in (|6.5| )) plugging (| 
into ( |6^ ) yields a quadratic equation for qg. Solving this equation we find 



Im — 3q 

Qs= \q 



m — 3q 



(6.10) 



The squared radius of the black hole is given by R^ = 2qh. At leading order the entropy 
and temperature of the mixture are given by 



,3 _l n l^-^1 



T 



s = 7rR^ = n\ 2\^ 

(6.11) 
\/2 fm-2,q\^ 
~K V 3 

As the m — 3q for the mixture is of order q'^ , it is convenient to define a shifted and rescaled 
mass variable 

q'^p = - (m - 3q) 
in terms of which 



3 3 

s = 7rq2 pi 

pi (6-12) 



We will see below that these results correctly reproduce the leading order thermody- 
namics of hairy black holes. 



6.4 Hairy Black Hole 

Once we have our solutions for hairy black holes from Appendix ^, the evaluation of their 
thermodynamic charges and potentials is a straight forward exercise. At low orders in the 
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perturbative expansion we find 
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(6.13) 
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+ e^ 



^ + OiR') 
6 



+ e^ [0(i?^)] + e^ [0(i?2)] + 0(e«) 
47rr = [4i? - 2R^ + ©(i?^)] 



+e^ 
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+ e'^ 



f + 0(fl^) 



+ e» [0(ii)] + Oie") 



These quantities obey the second law 

dM -3fidQ- TdS = 

where S = Entropy = ^R^ = ttN^R^. 

Upon setting e = in ( |6.13| ) we find we find a formula for the instability curve of the 
RNAdS black hole. Eliminating R we find that this curve lies along the curve 



m = 3q + 3q'^ + 0{q'^ 



(6.14) 



i.e. along the curve 5^ = 6 in the notation of subsection ^^. Note that, along this curve, 
H deviates from the value unity only at 0{R'^). 

We will now compute the entropy of the hairy black hole as a function of its mass and 
charge. Let us define a rescaled energy above BPS 



4 
3 



pq = - [m — 3q) 
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It may be verified that 



R' 



qVp[l + l{-2 + ^) + 0{q') 



Vp) 



{2 + ^-p)+0{q') 



(6.15) 



The temperature and chemical potential of the black hole are are given by 

1 
T = ^^[4 + q{6-5^) + 0{q^)] 

l, = l+q^ + Oiq') 



.16) 



while its entropy is given by 
S = ttR^ 

+ 0(q^ 



3 3 

2 7rp4 



i^3g. _ ^, , 100-76^+13p 
l + -i-2 + ^p)+q 



.17) 



Note that ( 6.16| ) and (|6.17| ) agree with ( |6.12 ) at leading order. It follows that, atleast 
for thermodynamical purposes, the hairy black hole solution may be regarded as a non 
interacting mix of a small RNAdS black hole and the soliton, at leading order. 



7. Hairy Rotating Black Holes 

7.1 Thermodynamics of small Kerr RNAdS Black Holes 

Large classes of explicit Kerr RNAdS black hole solutions have been presented (and their 
thermodynamics worked out) in ||l|, |2|, ^, ^, ^, ^, 0]. In this subsection we focus on the 
special case of small near extremal black hole with self dual angular momenta and all three 
charges equal. Concretely, we study black holes with mass M = N'^m, charge Q = N'^q 
and angular momentum J = N'^q^j ^^. We will take q to be small, but allow j to be 
arbitrary. All our formulae below are presented in a power series expansion in q but are 
exact in j. 

The extremality curve for the black holes we study is given by 

me.t{q,j) = 3^ + (3 + ^) g' + (-6 + ^/ - ^A q^ + 0{q^) (7.1) 



In contrast the BPS bound for the theory is 

rriBPsiqJ) = M + 2q^j 



Note that 



rriext - ruBPS = q 



: 0- - 3)^ 



+ 0(g3 



(7.2) 



(7.3) 



■^^ Concretely, J is the value of Jz in one of the two SU{2) factors of 50(4). As in the earher part of this 
paper, the mass M is normahzed to agree with the scahng dimension of dual operators, while the charge Q 
is normalized to be unity for a complex chiral field Z. 
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In particular the mass of the extremal black hole always > that of the BPS black hole. 
The extremal black hole is also BPS only when 

j = 3 - 2g + 3g2 + 0{q^) (7.4) 

As in previous sections, we will be interested in black holes whose energy deviates 
from extremality only at order 0{q^). To focus in on these energies, it is convenient, as in 
previous sections, to define a shifted and rescaled energy variable 5"^ by 
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As above we work in a power series expansion in q but our formulae are all exact in 6. 
With this notation, the thermodynamical potentials of small near extremal RNAdS black 
holes is given by 
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(7.6) 



2V-2j^^^±f^)q^Oiq^) 



We are particularly interested in two different two parameter surfaces in this 3 dimen- 
sional space of black holes. The first of these surfaces is the subspace of extremal black 
holes given by ([T^). Specializing to this surface we find 



(7.7) 



The second surface of interest, as in previous sections, is the subspace of black holes 
with ^ = \. We may solve for 6 in terms of q in order to set /i = 1; the solution is given by 
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Substituting this value of 5 in 



n = | + (^(?±^,. 



486^/6(j2 - 9) 
we find 
2^^ /i(-18 + j2) 



q + 0{q^ 



(7.8) 
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Notice that the /x = 1 surface intersects the extremaUty surface along the hne of 
super symmetric black holes. Indeed, plugging the relation (7.4) into ( |7.7D we find 






(7.10) 



In fact it may be shown that the equation /x 
holes ^"^ 



ri = 1 is exact for super symmetric black 



7.2 Hairy Rotating Black Holes as a non interacting mix 

We have seen in the previous section that the thermodynamics of Hairy AdS black holes 
is very simply reproduced, at leading order, by a simple physical picture of the black hole 
as a non interacting mix of the RNAdS black hole and the soliton. In this section we will 
simply assume the same is true of charged rotating hairy black holes. In other words we 
assume that there exist charged rotating black holes whose thermodynamics, at leading 
order, is reproduced by an arbitrarily weakly interacting mix of the spherically symmetric 
soliton of previous sections ^^ , and the small spinning RNAdS Kerr black hole. As in the 
previous section, the condition for thermodynamical equilibrium of this mix is simply the 
requirement that the Kerr RNAdS black hole, that participates in this mix, has n = 1. 
The charge, angular momentum and energy of such an equilibrated mix is given by 



q = Qb + Qs 



jq^ = ml 
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3qb + -{jl + 9)ql + 0{ql) 



(7.11) 



+ ^qs 



The three equations ( [7.11 ) may be used to solve for qi,, qs and jb in terms of q, j and 
At leading order we find 
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(7.12) 



^^We thank Seok Kim for explaining this to us. 

^*It should also be possible to study rotating black holes in equilibrium with solitons made out of the 
condensates of other supergravity modes, e.g. modes of the graviton. At small total angular momentum, 
however, the only modes of the chiral ring satisfy the thermodynamical requirement of non interacting 
equilibrium. This changes at angular momenta exceeding unity at which point equilibrium between spinning 
black holes and the solitons of this paper becomes impossible (see below), but rotating black holes could 
presumably equilibrate with other solitons. We thank H. Reall for a discussion on this point, and leave the 
discussion of black holes immersed in other solitonic backgrounds to future work. 
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It then follows that the entropy, temperature and angular chemical potential of the 
mix, are given at leading order by 
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Note that the hairy black holes of this subsection exist only for j < 3 (this follows 
because j^ > j whenever g^ > (as we have assumed), but black holes with j^ > 3 have 
negative temperature and so a naked singularity). When this condition is satisfied, they 
exist provided 



3q + 2jq^ <m<3q+-{f + 



^jHf+9)q' 



(7.14) 
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The range of this existence interval shrinks to zero as j approaches 3. 

At the upper end of this mass range we find, from ( 7.12| ), that qi, - 
qs = 0; i.e. the hairy black hole reduces to a pure RNAdS black hole. At the lower end 
of this range (i.e. when the system charges satisfy the BPS bound m — 3q = 2jq^) on the 
other hand ( |7.15| ) reduces to 
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(7.15) 



q\l 



At this lower end, then, the temperature of the hairy black vanishes. At this end the 
hairy black hole is neither pure soliton (unless j = 0) or pure black hole (unless j = 3; 
at which value the upper and lower end of (|7.14| ) coincide) but a mix of soliton and black 
hole. The special feature of this mix is that jh = 3, so that the participating black hole 
is super symmetric. At the lower edge, then, the system is a weakly interacting mix of a 
supersymmetric black hole and the supersymmetric condensate, and is itself supersymmet- 
ric. This is of course intuitive. An extremal hairy black hole is given by a non interacting 
mix of an extremal Kerr RNAdS black hole with ^ = 1 and a soliton. But, as we have seen 
above, extremal Kerr RNAdS black holes are BPS at /i = 1. It follows that (within the 
approximations of this subsection), that the extremal hairy black hole is also BPS. 

Note that the approximations of this subsection predict a two parameter family of BPS 
hairy black holes (the entire lower edge of ( 7.14D ) in comparison with the one parameter 
set of BPS Kerr RNAdS black holes. These black holes constitute the lower end of the 
mass range ( [7.14] ). 



54 



To end this section we should emphasize that ah the formulae of this section are 
predicated on the assumption that hairy black hole thermodynamics may be reproduced, 
at leading order, by a non interacting mixture of a Kerr RNAdS black hole and a scalar 
condensate. While we have checked by explicit calculation that this guess is true in the 
absence of rotation, we have not yet performed any such check for rotating black holes. 
The reader should, consequently, regard the formulae of this section as conjectural. As 
we have remarked in the introduction, it may well prove technically possible to verify (or 
diversify) the predictions of this subsection by an explicit perturbative calculation of the 
sort presented earlier in this paper. We leave this effort for future work. 



8. Discussion 

In this paper we have studied charged black holes in global AdS spaces. We have focused 
mainly on spherically symmetric black holes with equal diagonal S0{6) charges. At small 
values of the charge we have demonstrated that the spectrum of AdS hairy black holes 
extends all the way down to the BPS bound. We have also conjectured that this result 
applies at all values of the charge. The evidence for this last conjecture is not yet over- 
whelming, it would certainly be worthwhile to gather other evidence (e.g. from numerical 
solutions of the relevant differential equations) to support or refute this claim, as also to 
verify the precise nature of the large charge singular super symmetric solutions proposed in 



subsection 5.3.2 



The special singular solution S, discussed in section |5|, appears to play a very special 
role in the space of supersymmetric solutions of the theory. It would be interesting to 
study the near singularity geometry of this solution in more detail; in particular it would 
be very interesting to determine the 10 dimensional lift of this solution, as well those of the 
a = 1 singular solutions of section |^. It is conceivable that these singular five dimensional 
solutions have a regular ten dimensional lift. If this is not the case the ten dimensional 
perspective should at least yield valuable insight into the nature of the singularities in these 
solutions. It may also be worthwhile to re investigate the nature of a = 2 supersymmetric 
solutions from the ten dimensional viewpoint. 

Let us momentarily turn to the consideration of black branes in asymptotically Poincare 
AdS spaces. These solutions may be obtained from the large charge limit of the black 
hole solutions studied in this paper. RNAdS black branes solutions exist at all values of 

4 

Pe > CePq with Ce = f VT^'^. These solutions are presumably unstable to the condensation 

4 

of the (j) field for energy densities pe < CgPq for some Cg > Cg. If hairy black branes at 
a given charge density have a lower bound on their energy density, the equation for this 

4 

lower bound must also take the form pe > CmPq with Cm < Cg. It is clearly of interest to 
know the value of Cm- Our conjecture that hairy black holes descend down to their BPS 
bound even at large charge amounts to the prediction that Cm = 0, i.e. there is no lower 
bound for the mass density of hairy black branes, at any fixed finite charge density. This 
phenomenon is not without known precedent; as we explain in appendix^, a similar result 
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is true of one charge (rather than 3 equal charge, as studied in this paper) RNAdS black 
branes. 

On a related note, the recent paper |3(:] has presented an analytic determination of 
the constant Cg (and the hairy black brane solution in the neighborhood of this critical 
density) in a system with some similarities to the one studied in this paper. In particular 
the scalar field in [^6| has the same value of the mass as in the current work, but carries 
infinite charge (this is the so called the probe approximation that ignores the backreaction 
of gauge dynamics on the metric) . It would be interesting to investigate whether the results 
of [^ can be generalized to the study of our Lagrangian (|l.2D . 

It would be of interest to extend the results of this paper to the case of charged black 
holes with three unequal SO{6) charges, and particularly to study special limits in which 
one or two of these charges are turned off. In particular, in the special case where only 
one SO (6) charge is nonzero, ordinary RNAdS black holes extend all the way down to the 
BPS bound, so there is no pressing reason for these solutions to exhibit a superradiant 
instability, or for one charge hairy black brane solutions to exist. This question deserves 
further investigation. 

It it would be fascinating to verify the correctness or otherwise of the tentative pre- 
dictions of section [^ for the spectrum of small rotating black holes. In particular, if new 
physically acceptable super symmetric hairy black holes really do exist, it would be fasci- 
nating to determine them and to analyse their properties. 

As we have mentioned in the introduction, it is possible that the black holes we have 
constructed suffer from a superradiant instability towards the emission of the gravity modes 
dual to TrX^ + TrV^ + TrZ^ . It should in principle be straightforward (though perhaps 
tedious in practice) to check whether or not this is the case, by computing the imaginary 



part of the relevant quasinormal modes, along the lines of Appendix A of |14]. 

In this paper we have focused only on a very particular kind of hairy black hole; one 
whose condensate is the zero mode of a scalar field in AdS^ space. It is likely that there exist 
many different such hairy solutions with various different gravitational field condensates. 
In the current paper we have likely illuminated only a very small part of an intricate and 
fascinating structure of hairy black hole solutions in AdS'^ x S^. It would be fascinating 
to more fully uncover the structure of new black hole solutions in AdS^ x S^, and perhaps 
most importantly, to understand their properties directly from the dual Yang Mills point 
of view. 
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A. Results for the Perturbative Expansion of Hairy Black Holes 



A.l Far Field Solution 



Scalar Field: 
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Metric and Gauge Field: 
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1152r2(l + r2)^ 



48 (1 + r 



2^3 



log 



(A.2) 



A. 2 Intermediate Solution 



Co)(y)=l 

4>Tii)^y) =-V^ - 41og(i?) - 2 log (y2 - 1) - 2 

</)[^^4^)(y) =81og2(i2) + Uy^ + 81og (y^ - l) + ^^ + 26 J log(ii) + ^ _ i 

log (7/2 - 1) (2/ + 11^2 ^ 2 (^2 _ 1) log (^2 _ 1) _ ^ ^ ipg(4) _ ig 

y2-l 



4 _\ . .„. / + 5y^ + 4y2 



+ 
1 






-3y2 - 121og(i?) - 6 log (y2 - l) - 5) 



imid 



(y) 



55 



(5,0) Vi/; "2304 

w / ^ _ -813j/2 - 32521og(i?) - 1626 log {y^ - l) - 1091 

{5,2){y> - 6912 



(A.3) 
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,m^d,..^ (r-1)^ 



Metric and Gauge Field: 
/(To) (2^)= 4y4 

Co)(y)=o 

,^id. , _ {y^ - 1) [/ + (%' - 4) logjR) + (%^ - 6) log {y^ - 1) + log(4) - l] 
H2,2){y) - 4^4 

mid , ^ _y' [y' - 41og(i?) + 21og (y2 _ l) _ log(4) - 5] 



C2)(y) 



4(y2_l 



mid ,..^ _/ + (8y' - 4) log{R) + (4y2 _ 6) log (y^ _ l) + log(4) - 1 



8y' (A.4) 

Co)(y)=o 

j(4,2){y) - yq^ 

mid ( \ 



fTm^y)=o 

jmid 



mid , ^ _/ W - 41og(^) + 21og (y2 - 1) _ log(4) - 7) 



16(y2-l)^ 

.^rf , , _ / - 2j/^ + (8;/^ - 4) log(i?) + 2 (2y^ - 3) log [y" - l) + log(4) + 5 
^(4,2) 12/) - ^4^2 



A.3 Near Field Solution 

Scalar field 



C2)(^) = - 81og(^) - I log(^ + l)(log(^ + 1) + 4) - Li2(-z) - 21og(2) - ^ - 3 



,.„ . . 2(4 + 31og(2))z + TT^jz + 1) + 61og(2) + 12 
'^(3,2)W= ^(T^Y^ (A.5) 

+ |(log(z) + i^ - 2) log(z + 1) - 3 log(i?) + ^Li2(z + 1) 
o o 

(3z(2 + l) + 2)log^(z + l) 

Wz{z + 1) 

55 
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Metric and Gauge Field: 



.2 1 



fll4z)=-z'--{2z + l)\og{z + l) 

in I ^ _ (2Z + I)l0g(z + 1)-4Z 
^(2,-2)^^^- 32z2(^ + l)2 

A'C2,2)i.^)=\{-Z-\0g{z + l)) 

,„ _ z(4z + 2)+lo g(z + l)((z + 1) log(z + 1) - 1) (A.6) 

7(4,4) l^i — 



9it-2)(^) 



16(z + l) 

(5 - 12z)z2 + {A{z - 2)z - 3) \og{z + l)z + (3z(z + 1) + 1) \og^{z + 1) 
" 256z3(z + l)3 

z(z + 21og(z + l)-3) 



^(^'2)(^)= 64(z + l) 

B. Supersymmetric solitons in A(iS'5 x S^ and the planar limit 

B.l One-charge solitons 

In this paper we studied supersymmetric configurations in AdS^ x 5^ by considering a 
truncation of AA = 8 supergravity where all 3 charges corresponding to the three orthogonal 
planes of SO{&} are taken to be equal. The spherically symmetric supersymmetric solutions 
of this truncation are 1/8 BPS from the ten dimensional point of view. It is instructive 
to briefly review solutions of A/" = 8 supergravity where we turn on only one of these 
charges. Spherically symmetric supersymmetric configurations of this type are 1/2 BPS 



and belong to the LLM family |37]. We look for static, spherically symmetric 1/2 BPS 
solitonic solutions with the topology of AdS^ x S^ (so that in the LLM language they have 
the topology of a disk). These can be studied in a truncation of AA = 8 to f/(l)^ coupled 
to 3 hypermultiplets, by turning on only one of the U{\) charges and can be expressed in 
the language of |15, ^] as'^^ 



Ai = A2= 0, ^3 = h-^dt, Xx = X2 = /^^/^ Xg = /i^^/a (B.l) 



»2 



0, </.3 = 2 V(/i + p/iV2)2 - 1 



where Ai are the 3 t/(l) gauge fields, Xi the scalars in the vector multiplets constrained 
to satisfy X1X2X3 = 1 and </>i scalars in the hypermultiplets. This is a solution if h{^p) 
satisfies the equation 

(1 + p^K){?>U + ph!') = p[4 - (2/i + ph!f\ (B.2) 

Notice the different powers of h in the metric ( p.lj ) relative to the 3-equal-charge case (p.lj). 
Also, unlike equation (5.7), now ([B.2|) can be solved analytically and we find the following 



^^We use somewhat different notation with rthere = Phere, (-H'l)there = {H2)theri! = 1, {Hz)there = 



p'l.h.ere • 
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set of smooth solutions |15, |lq 



The parameter q corresponds to the charge. In contrast to the 3-charge case which was 
studied in the rest of the paper, now q is unbounded from above i.e. there is no maximum 
charge of smooth spherically symmetric 1/2 BPS solitons similar to ( [5.22 ). 



In particular we can try to recover a spacetime which is asymptotically the Poincare 
Patch by taking the large charge limit. We consider the following scaling 



t = k T, p = ku, Qi = k Xi, q = qQk 
T ,u ,Xi ,qo = const , fc — )• oo 

Then we find the following solution 



(B.3) 



, .. , .. ,. fl + l^V .. 



ds^ = -u^ ll + -f ] dr^ + —du^ + uM 1 + -^ I dxt 

A, =A2 = A^ = 



-^(1 + lf 



4go 



(B.4) 



4'3 



t\/4go + 



u^ 



which is asymptotically the Poincare Patch and has scalar fields turned on in the interior. 
Let us finally consider the non-extremal one-charge RNAdS black branes. We start 
with a non-extremal one-charge RNAdS black hole in global AdS. We have the following 
solution 

ds'^ = -H-'^l^fdt^ + H^'^{f-^dp'^ + p'^dnj) 



Ai=A2 = 0, ^3 = 2 . r^ "^ 



P^ + b 
X^ = X^ = H^'\ X' = F-2/3 (B.5) 

<Al = 02 = </'3 = 



whose charge is g = \/h{p + h) and mass m = |p + g. The black hole is extremal (and 
singular, i.e. the superstar) when /i = 0. We now take the scaling 



t = k ^ T, p = ku, ^i = k ^ Xi, b = b'k'^, fi = fi'k^ 
T ,u,Xi, b' , p = const , /c — )• oo 
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(B.6) 



Then we find the fohowing non-extremal charged brane solution 

ds^ = - (l + ^) "'' {-' + ^' - ^) dr^ + (l + ^) ''' [{-' + ^' - ^) " d-' + ^^-^^ 
Ai = A2 = 0, A3 = i ,^,,o dt 

/ Ui \ 1/3 / y \ -2/3 

X^=X2=h+— j , X3=h+— j , </>! = ,/,2 = 03 = 

(B.7) 

The energy density and charge density of this solution is (up to numerical factors of order 
one) 

Pe w At', Pq ~ \Ay 

We notice that 

pq 

we notice that by taking /x' small enough we can make this ratio as small as we like. Hence 
for any finite charge density, the energy density of one-charge RNAdS black branes can 
be made arbitrarily small i.e. the constant Ce mentioned in section ^ is Cg = for the 
one-charge case. 

B.2 Three-charge solitons 

A similar scaling can be performed for the 1/8 BPS (3-equal charge) soliton. For this one 
takes the following scaling of solution of the - form 

t = k~ T, p = ku, fli = k~ Xi, a = 2JTink 
T ,u ,Xi ,qo = const , A; — )• 00 

For large a the solution ( |5.37 ) becomes a good approximation and from ( |5.1[ ) we end up 



with the following hairy, asymptotically Poincare Patch solution of (2.4) 



A = , 6 '^° 



ly/4qo + 



u^ 



It would be interesting to explore whether the solution ( |B.9| ) has any relation to hairy black 
branes near the BPS limit. 

C. Some numerical results 
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